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Abstract 


The first phase of the study of the performance of a wormgear 
transmission is reported. In this phase the work included the selection of a 
double -enveloping wormgear type, and its dimensions, suitable for use in 
helicopter transmissions; the 3-D graphics representation of the selected 
wormgear using the I-DEAS software; the analysis of the kinematics of meshing; 
the analysis of load sharing among the meshing teeth; and the implementation 
of the analyses in a computer program. The report describes the analyses, 
their results , and the use of the computer programs . 
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1 Background and Objectives of the Study 


The wormgear is attractive for use in helicopter transmissions because 
it can accomplish two functions in one stage, viz. transmitting torque from 
the horizontal engine shaft to the vertical rotor shaft, and simultaneously 
realizing the required speed ratio. Thus, the use of the wormgear can result 
in compact transmission design with significant weight saving. However, the 
wormgear has a serious drawback. Because worm and gear teeth are in pure 
sliding contact, the sliding friction causes excessive power loss and renders 
the efficiency of power transmission low. Recently, Chaiko (1990) advanced 
the concept of using hydrostatic lubrication to remove the sliding friction. 

He considered a wormgear transmission of realistic dimensions and under 
realistic operating condition, and found that the efficiency of the 
transmission could reach 99.14%, and its weight could be only 41% that of the 
baseline transmission. These figures convincingly demonstrate that the 
hydrostatically lubricated wormgear can be a candidate of the advanced 
helicopter transmissions. To further investigate the performance of this type 
of transmission, an analytical tool capable of generating quantitative data 
for a wide range of operating conditions is needed. Hence a research proposal 
was submitted, and a grant was subsequently awarded (NAG3-1316), to develop 
such a tool. 

The project consists of three years of study with the following 
obj ectives : 

1. Year 1 is to be devoted to the selection of a double -enveloping wormgear 
type suitable for use in helicopter transmissions, the analysis of the 
kinematics of meshing, and the analysis of load sharing among the 
meshing teeth. 

2. Year 2 is to be devoted to the selection of a hydrostatic lubrication 
configuration on the gear surface, and a comprehensive analysis of the 
hybrid (hydrostatic and hydrodynamic) lubrication at the contact based 
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on the surface geometry and running conditions obtained in Year 1 . 

3 . Year 3 is to be devoted to the analyses of accessory bearings in the 

wormgear transmission, and the assembly of the obtained information in a 
computer program package for the design of this type of transmission. 

The present report describes the work done in Year 1. 


2 Wormgear Geometry 
2.1 Selection of Wormgear Type 

Wormgears are used to transmit power between two non-parallel non- 
intersecting shafts. The angle between the shafts is usually 90 deg. There 
are two common types of wormgears: single -enveloping and double -enveloping. 
Fig. 2.1 shows a single -enveloping wormgear where the gear envelops the worm. 
The double -enveloping wormgear is shown in Fig. 2.2, which is like the single- 
enveloping wormgear except that the worm also envelops the gear. 

The double -enveloping wormgear has two contact lines on the mating 
surface during most of the meshing period. In comparison, the single- 
enveloping wormgear has only one contact line. Hence, the double -enveloping 
wormgear is expected to have superior contact properties. The usual way to 
manufacture a wormgear is first milling a worm, then cutting the gear with a 
gear hob which has the same shape as the worm. The worm and gear hob are 
usually ground for high precision applications. The grinding wheel can be a 
disc (Shen, et al., 1983), Fig. 2.3, or a cone (Simon, 1989), Fig. 2.4. The 
worm grinding process is called the first enveloping process; and the gear 
cutting process is called the second enveloping process. 

In the first enveloping process the grinding wheel actually rotates 
around a base circle, with angular velocity V> 2 > while the worm rotates about 
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its own axis with angu.l 3 . 1 r velocity '0^. The distance between the won axis and 
the base circle axis is denoted as a 0 , and the ratio ij > x to t/> 2 is denoted as i 0 . 
In the case of the disc grinding wheel, the grinding plane may be parallel to 
the base circle axis or inclined to it at an angle {}. The center distance a 
and speed ratio i of the wormgear are actually the respective quantities used 
in the second enveloping process. If a=a 0 and i=i 0 , the wormgear is said to 
be of the conventional form; otherwise it is said to be of the modified form. 


In helicopter applications the torque load and the speed ratio are high. 
To accommodate such stringent operating conditions the double -enveloping 
wormgear was selected. This is not only because the double -enveloping 
wormgear has a larger contact region but also because its contact pattern is 
favorable to the establishment of fluid film lubrication. Further, for ease 
of manufacturing the disc grinding method was chosen. The analysis presented 
below is based on the conventional form, but it can be easily extended to the 
modified form. 

2.2 Layout of Wormgear Dimensions 

Wormgear dimensions are usually determined by following the gear 
standards (e.g. AGMA, 1965), which were developed based on many years of 
experience. The procedure typically starts with a known speed ratio and a 
known center distance. For helicopter applications, the suggested 
transmission power output is 4000 kW , the speed ratio is 110, and the center 
distance is about 1500 mm . ( Chaiko , 1990) These parameters are beyond the 
range covered by the AGMA standards. Hence, while the standards were followed 
as closely as possible, some modifications had to be made. Table 1 lists the 
given and the calculated dimensions and parameters, along with the formulas 
used for the calculation. 

Common terminologies used to describe the double -enveloping wormgear are 


3 



shown in Fig. 2.5. Another term that should be introduced is the mid-plane, 
which refers to the plane passing through the worm axis and perpendicular to 
the gear axis. The pitch circle of the gear is measured in the mid-plane. 

The pitch circle of the worm is measured at the center of the worm. In the 
double -enveloping wormgear the worm circular pitch and the gear circular pitch 
are equal. The base circle is used to define the worm profile in the mid- 
plane. The normal pressure angle is always chosen to be 20 deg. 

Several of the selected quantities in Table 1 should be explained. A 
right hand worm was selected. The pitch diameter of the worm d^^ was chosen to 
be 300 mm because of strength and rigidity considerations. The number of 
meshing teeth N p was chosen to be 10 instead of the usual 4 or 5 in order to 
reduce the load level on the meshing teeth. The choice of (3 = 5 deg was a 
compromise between acquiring a better contact condition and maintaining an 
acceptable tooth thickness at the top of the worm thread. 

It should be noted that for the double -enveloping wormgear a different 
choice of dimensions can result in very different surface shape and contact 
pattern. The results shown below are based on the data of Table 1. 

To better visualize the selected wormgear, a software named I-DEAS was 
used to generate its 3-D graphics, Figs. 2.6 to 2.9. Fig. 2.6 shows the 
engaged wormgear. The worm and gear are shown separately in Figs. 2.7 and 
2.8. Fig. 2.9 depicts a close-up view of the gear teeth. 

Details of the I-DEAS program that generated these graphics, and the 
procedure to run the program, will be described in Section 5.1. 
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3 Contact analysis 


3.1 Method of Enveloping 

The geometric and kinematic properties of the wormgear were analyzed 
with the method of enveloping. To describe the method, let us consider two 
rotating machine components, m and k, which are in meshing. Let us attach a 
coordinate system S m (x m ,y m , z m ) to the machine component m; and likewise attach 
a coordinate system S k (x k ,y k ,z k ) to the machine component k. The two machine 
components are rotating with the angular coordinates and with respect to 
some stationary reference frame. Since the two components are in meshing, 
there must be some definite relationship between V m and V»k* Hence only one of 
them is independent, and either can serve as the running parameter. Let tJ> 
denote this running parameter. The surface of machine component m can be 
expressed as f (x^y^ z m )-0 in S m . By coordinate transformation, f=0 becomes a 
family of surfaces, F(x k ,y k ,z k ,^)=0 in S k . Then the surface of the machine 
component k is the envelope of F=0, viz. 

F{x k ,y k ,z k . tJO=0 (1) 

dF/d\ |r=0 

If V’k/V’m = constant, the surfaces of the two meshing components are said to be 
conjugate surfaces (Litvin, 1989) . 

While the idea was simple, the algebraic manipulations needed to carry 
it out were quite clumsy. Recently, Shen, et al. (1983) reformulated the idea 
using a set of operators. Because many of the operators involving the 
transformation of coordinates can be done once for all, the method has become 
simple as well as elegant. In the special case of conjugate surfaces, the 
method is even more convenient to use. This will be demonstrated below. The 
reader is referred to Shen, et al. (1983) for the details of the method. 
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3.2 Coordinate Systems and Operators 


To carry out the analysis, four right hand coordinate systems, S lf S 2 , 

Sj and S n> are set up. See Fig. 3.1. 

Sj and S n are stationary systems. Sj(Xj,yj,Zj) is located at the worm 
position. Xj is along the center line of the wormgear and points to the gear, 
and Zj axis is along the worm axis. S^x^y^Zn) is located at the gear 
position. Xjj is along the center line but opposite to Xj, and z n is along the 
gear axis . 

S x and S 2 are moving systems. S 1 (x 1 ,y 1 ,z 1 ) rotates with the worm and its 
z^ axis always coincides with the Zj axis . The worm rotation angle is the 
angle between x x and Xj . S 2 (x 2 ,y 2 ,z 2 ) rotates with the gear and its z 2 axis 

always coincides with the z n axis. (In the first enveloping process S 2 turns 
around the base circle with the grinding wheel.) The rotation angle ip 2 is the 
angle between x 2 and x^. Thus, S x is fixed to the worm and S 2 is fixed to the 
gear; and the worm and gear rotate about Sj and S n , respectively. Initially 
S x coincides with Sj, and S 2 coincides with S n . 

Let A^ (km) be the coordinate transformation matrix from S m to S k , where 
the subscript ( indicates that ip is the running parameter of the coordinate 
transformation. If the running parameter is understood, the subscript is 
dropped to simplify the notation. In addition, if the transformation is 
between two stationary systems, then no running parameter is involved and 
there will be no subscript. The coordinates in S m can be transformed into S k 
as : 


x k 



y k 

= Ai, M 

y m 

Zk 


Z m 

i 


1 . 


( 2 ) 
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and a function in S ro can be transformed into S k as: 

A im f(x„.y m .z„) =f <3) 

The coordinate transformation matrices from Sj to S x , from S 2 to S n , and 
from S n to Sj are given as: 

cosily sini^ 0 0 

A iW = - sin ^i 0 0 (4) 

0 0 10 

0 0 0 1 , 

cost|r 2 -sini|r 2 0 0' 

A (^) = sin *2 cosy 2 0 0 (5) 

0 0 10 

0 0 0 1. 

-10 0 a 0 

0 0-10 (6) 

0-100 
0 0 0 1 . 

The transformation matrix from S 2 to S x may be written as A^ 12J =A (1 ^ ) A ( ^ n) A (n2) , 
or : 

-cosijrjcosilr.j cosies intjr 2 -sini^ a^costl^ 
sin^cos^ -sin^sin^ -cosi]^ -aj^sintlr., 

-sini|r 2 -cosi|r 2 0 0 

0 0 0 1 

Note that the relation A <njk) =(A (km) ) -1 holds. 

In addition to the coordinate transformation matrix, two other operators 
are defined. One is the differentiation operator: 





7 



( 8 ) 



Dj=E> f L> t = 


a 2 

ch|r 2 


The other is the similar differentiation operator: 

^ 1} = Ai 12) D^Ai 21) (9 

Z^ 2> = 

In the above formulas, the subscript ( )^ may be dropped if no confusion can 
arise . 


Several needed properties of these operators are given below: 

DA {km) f{x m ,y m ,z m ,\ Jr) (10) 

=A (km) Df+ ( A {km) f x ) DA (km) x m + (A {km) f y ) DA {km) y a + U {ka) f z ) DA {km) z m 


DA {ij) f=A ii i ) DV ) f 
A li i ) Df=D U) A li 3 ) f 


(ID 




( 12 ) 


where 


r _ 3/ f _ df f - 3f 
*■ dx m ' z y- dy m ' 2 - 3z„ 


(13) 


In the above formulas, k and m can be 1, 2, j and n; but i and j can only be 1 
and 2 . 


To distinguish between the first and the second enveloping processes, tJj 
is used as the running parameter for the first enveloping process and 6 is 
used for the second enveloping process . The order of operations with respect 
to different running parameters may be interchanged, e.g. 
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( 14 ) 


The derivation of the above and other properties of the operators can be found 
in Shen, et al. (1983). 

It should be noted that different symbols, i 0 and i, are used to denote 
the respective speed ratios in the first and the second enveloping processes, 
i . e . 


. = = 

■ i ° t * 

*2 ^2 

(15) 

02 02 

As said before, if i 0 - i and a 0 =a the wormgear is of the conventional form; 
otherwise it is of the modified form. 

3 . 3 Worm Surface 

The worm surface is formed by a grinding wheel in the first enveloping 
process . The surface of the grinding wheel is called the original generating 
surface. For a disc grinding wheel this original generating surface is a 
plane. For a cone grinding wheel the original generating surface is a line on 
the cone. 


The original generating surface can be expressed in the S 2 system as 
f (x 2 ,y 2 ,z 2 )=0. In the first enveloping process S 2 rotates about z 2 with the 
rotation angle 0 2 , &Tid rotates about z^ with the rotation angle 0^. As a 
result of rotation, the original generating surface forms a one -parameter 
family in S 1} F(x 1 ,y 1 ,z 1 ,0)=A^ (12) f (x 2 ,y 2 ,z 2 )-0. Because 0 x /0 2 = constant, the 
running parameter 0 can be either 0^ or 0 2 . Then, the envelope of the family 
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of surfaces is the worm surface, which is given as F=0 and D^F=0. 

Considering only the case of the selected wormgear, the original 
generating surface is a plane in S 2 , which is: 

f (x 2 ,y 2 ,z 2 ) =y 2 -z 2 tanp+r b =0 (16) 

After coordinate transformation into S-^ the one-parameter family of planes 
becomes : 

F=A^ Z) f=x x (cosT|r 1 sinilr 2 +tanPsimlr 1 ) 

+y 1 (-sinTlTj^sin^+tanpcos^j) ( 17 ) 

-j^cos^- (a 0 sim|r 2 -r h ) =0 

Then, the worm surface may be represented by the following two equations: 

F = Xi ( 003^3^1^2+ tanps inti) + ( -sint a sint2 + tanPcosti> 

-^cosil^-UoSin^-r*) = 0 

D^F = x 1 (cosTlriCOSt2 +i otanpcosti“i 0 sin ^i sini l r 2^ (18) 

+ y 1 {-sin^ 1 cos^ 2 -i 0 tan^sin^ 1 ~i 0 cos^ x sim^ 2 ) +z 1 sini|r 2 -a 0 cost|r 2 =0 


Using the relation Vh/V^lo and eliminating Vh and V *2 from Ec l s - 18 » an explicit 
equation, which is nonlinear, of the worm surface can be obtained. 

For a given value of V>i» Eqs. 18 represent the equation of a straight 
line. Recall that the first equation is the equation of the original 
generating plane expressed at the position ti* Hence the straight line is on 
both the generating plane and the enveloping surface. This straight line is 
the contact line between the worm and the disc grinding wheel at the position 
jp x . It can also be seen that the worm surface is the set of all the straight 
contact lines. Fig. 3.2 depicts this situation. 
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According to differential geometry, if the planes in a one -parameter 
family are not parallel and there is not a point common to all the planes, 
then the envelope of this family is the tangent surface of its edge curve 
(Kreyszig, 1959). In Fig. 3.2 the curve labelled T is the edge curve. The 
edge curve is given by these equations : 


F = x 1 (costy 1 sini|r 2 +tanPsini|r 1 +y 1 ^-sintlr 1 sintlr 2 +tanpcosi|i 1 -.z 1 cosTlr 2 -(a 0 *sini|r 2 -r Jb ) = 0 
D^F=x 1 (cosi|r 1 cosiJr 2 +i 0 tanPcosT|r 1 -i 0 sin4r 1 sin^ 2 ) 

+ yit-simlrjcos^-iotanpsinijji-iocos^sin^) +z 1 sim|r 2 -a 0 cosi|r 2 =0 
DfD^F= x 1 (iocostjj 1 sini|r 2 +iotanpsini|r 1 +2*i 0 sim|; 1 cosiir 2 +cosilr 1 sim|r 2 ) ( 19 ) 

+ y 1 (-iosim|r 1 sinijr 2 +iotanpcost 1 + 2 *i 0 cosi|r 1 cosilr 2 -sini}r 1 sini|r 2 ) -z 1 cost 2 -ao sin ^ 2 = 0 
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The edge curve divides the envelope surface into two branches , of which 
one is the worm surface. If the worm thread domain determined by the selected 
dimensions intersects the edge curve, the surface of the disc grinding wheel 
would extend over the edge curve. The crossed-over portion of the grinding 
surface would undercut the worm thread. Undercutting the thread would 
decrease the length of contact lines, decrease the root thickness of the worm 
thread, and consequently decrease the load carrying capacity of the thread. 
Hence, wormgear dimensions should be chosen properly so that the worm thread 
domain is away from the edge curve. 

3.4 Gear Surface 

The same procedure applies to the second enveloping process. In this 
process the gear hob has the same shape as the worm surface. Hence, the 
generating surface is given by Eqs. 18 in Sj_. As the result of rotation, this 
generating surface forms a one -parameter family of surfaces in S 2 . The 
running par am eter for this process can be either 9 X or 0 2 , while ip becomes a 
geometric parameter already in existence in the worm surface equation. The 
procedure to obtain the gear surface may be as follows: By eliminating ip from 
Eqs. 18, one obtains an explicit equation of the worm surface. This equation 
may be written as : 

F*(x 1 ,y ir z 1 )=0 (20) 

The one -parameter family of surfaces in S 2 is obtained by coordinate 
transformation : 

f <2> U 2 ,y 2 ,z 2 ,e)= a} 21> F’= 0 

The envelope of the family: 

f <2) = 0 
D 9 f U) = 0 

is then the gear surface. 


( 21 ) 


( 22 ) 
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But it is not easy to obtain F*=0. Neither is it necessary. Assume 
that is obtained by solving the second equation of Eqs. 18, 

D^F=0. Substituting it into the first equation of Eqs. 18, F=0, yields the 
explicit worm surface equation 

F*= F {x x .y 1 ,z lt y(x 1 ,y x ,z 1 )) = 0 (23) 


Then, the equation D 0 f ( 2 ) =D 0 a/ 21 ) F*=O can be written as: 

DqAq 21) F*=Aq 21) D 6 F* 

+ US, 1211 F' Xi ) D^X 1+ (^ (211 FV,> ^ *V D »^ 21) Z l = 0 


(24) 


But from Eq. 23, 


F\ = 


(25) 


Then, since D^F-0 and Eq. 23 does not contain 6, Eq. 24 reduces to: 

DqAq 21 ^ F* = iAj 21> F Xi )D^ 21) X 1 +^ ) F yt )D^ 1 W 1 +{Ai Z1, F Xi )D l ,Aj™Z 1 =D l Ai 21 'F=0 (26) 

which, on multiplif ication by A fl (21) A 0 (12) , becomes A e (Z1) D e (1) F=0 . Thus, the gear 
surface is given by: 


Aq 21) F = 0 
Aq 21) D^F = 0 

Ae (21) £b {1) F = 0 


(27) 


which can be written as: 


a il a i2 a i3 a X4 
a 21 a 22 a 23 a 24 

a 31 a 32 a 33 a 34 
0 0 0 1 



*2 


0 ‘ 


y 2 


0 


Z 2 


0 


1 . 


.1 


( 28 ) 


where 


13 



a 11 =-sin4r 2 cos0 2 cos (tI^-Ox) -tanpcos0 2 sin (^i-Ox) +cos^ 2 sin0 2 
a 12 = sini|r 2 sin0 2 cos (tyx-Oz) +tanPsin0 2 sin(i|r 1 -0 1 ) +cosi|r 2 cos0 2 

a 13 = sin^ 2 sin(ilr 1 -0 1 ) -tanPcos (lifi-Oi) (29) 

a 14 = a [sinilr 2 cos (i|rx-6x) +tanPsin(i|r 1 -0 1 ) ] -aoSin^+rj, 


a 21 = sint|f 2 cos0 2 sin(ilri“0i) -tanPcos0 2 cos 
- [sini|f 2 sin0 2 +cosi|r 2 cos0 2 cos <i|rx - 0i) ] / i 0 
a 22 =-sint|r 2 sin0 2 sin (^-©x) +tanpsin0 2 cos (i|rx-0x) 

- [sinijr 2 cos0 2 -cosi(r 2 sin0 2 cos (t^-Ox) ] /i 0 
a 23 = sint|r 2 cos (i^-Ox) +tanPsin (^-©x) + [cosi}r 2 sin(i|rx-0x) 3 /i 0 
a 24 =-a [sint 2 sin(^x~0i) -tanpcos (t|Tx-0x) ~cosiJr 2 cos (\Jrx~0x) /i 0 3 -a 0 cos*|r 2 /i 0 


a 31 = a 12 -i* [sirulr 2 cos0 2 sin(ilrj-0i) -tanpcos0 2 cos (^x-0x) ] 
a 32 =-axx+i* [sini}r 2 sin0 2 sin (t|Ti-0x) -tanPsin0 2 cos (i|*x-0x) ] 
a 33 = [siniJi 2 cos (t|ix*~ ©x) ^tsnPsiii (31) 

a 34 = a*i* [sin\Jr 2 sin (t|Tx“0x) -tanpcos (^i~0x) 1 


In Eqs. 27 there are three equations and five unknowns. Using the 
relations ^ 1/^2 = *■<> and d i/ e z 3=1 the explicit form of the envelope can be 
obtained by eliminating i> and 9. Eqs. 27, along with the relation limiting 
the gear tooth domain, then provide the shape of the gear tooth. 

Consider the case of the modified form. For a given value of 9, Eqs. 27 
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produce a curve. The curve is the contact line between the worm and the gear 
at this 9 position. The curve goes over the gear surface twice as shown in 
Fig. 3.3. For fixed values of V and 9 , Eqs . 27 represent a point on the 
contact line . 

Consider the case of the conventional form. For a given 9 , if one 
chooses Tp=9 , the third and the second equations of Eqs. 27 are identical. 

Then Eqs. 27 reproduce a straight contact line of the first enveloping 
process. This kind of contact lines is called the first enveloping contact 
lines (the solid lines in Fig. 3.4). The set of all these straight lines 
reproduce a part of the original generating plane. For a given 9, let 9 
Eqs. 27 give a curved contact line, called the second enveloping contact line 
(the dashed lines in Fig. 3.4). The set of all the curved contact lines form 
the curved part of the gear surface. Thus, in the case of the conventional 
form, the gear surface consists of a planar region and a curved region. For a 
certain range of 9 values these two regions overlap each other. The 
overlapped region is shown in Fig. 3.4 where the dashed lines and the solid 
lines cross. This situation is impossible to occur in reality. From a view 
of the (y 2 ,z z ) plane it is clear that the curved region is farther away from 
the solid interior than the planar region. Hence, the curved region is 
fictitious, which is actually hobbed off during manufacturing. Analytically 
the principal relative curvature is used to sort out the fictitious part 
(Shen, et al, 1983). The relative curvature is the difference of curvatures 
of the worm and gear surfaces at the contact point. Because the gear tooth 
envelops the worm thread, at the contact point the curvature of the gear tooth 
must not be greater than that of the worm thread (absolute values). Hence, if 
the relative curvature is negative, the produced gear surface is not real. 

The derivation of the expression for the principal relative curvature is given 
in Appendix 3. 

From the above description it can be seen that, in the double -enveloping 
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wormgear, the worm and gear are in line contact and there are usually two 
contact lines on each pair of the meshing teeth. 


3 . 5 Tooth domains 

After the equations of the worm and gear surfaces are obtained, the 
tooth domains need be specified to determine what portions of these surfaces 
are the actual worm thread and gear tooth. Table 1 can be used to set up the 
following geometrical limitations on the worm and gear surfaces. The worm and 
gear tooth domains are shown in Fig. 3.5. Several new symbols that appear in 
the following equations are also defined in this figure. 

The maximum and minimum radii of the worm thread are: 

a- ^(<V 2+c)2 ~ z i 

a ~ V (c V 2+c)2 " z i 

r v = ^ 32 ' ) 

Z w Kln ^ r w MUt 


where r w is the radius of any point on the worm surface. 

The maximum and minimum radii of the gear tooth are: 
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R Gm =a- y /(d 1 /2) 2 -{b 2 /2p 


r^= min {(a+Aa) - (z 2 /2) 2 , R ^ } 

r^ta+Aa) - y/h^- (z z /2) 2 

r G = y/xZ+yi < 33) 

r <W r d* ^ 


where r G is the radius of any point on the gear surface, and 

h ^ a - ^d £ j2) 2 -{bj2) 2 

h mln = a - J{d a j2) 2 -{bj2) 2 (34) 

A a = (WV 2) " a 


3.6 Properties Needed for Lubrication Studies 

To carry out the lubrication analysis several geometric and kinematic 
properties are needed. These include the shape of the two contacting 
surfaces, the relative velocity between the two surfaces at a contact point, 
and the normal forces acting on the meshing teeth, etc. 

3.6.1 Surfaces in Stationary Systems 

Eqs . 18 and Eqs . 27 are obtained in the moving coordinate systems, where 
the worm and the gear appear stationary. To see the moving worm and gear 
surfaces, these equations must be expressed in the stationary systems. At a 
given instance , or for a given gear position 6 , there are two ways to obtain 
the worm and gear surfaces in a stationary coordinate system. One way is 
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calculating the coordinates in the moving system first, then transforming 
these points into the stationary system. The other way is transforming the 
equations into the stationary system first. In either way the coordinate 
transformation matrices are specialized to the given position. To distinguish 
this position from the running parameter 9 , it is denoted as 9 *. Thus the 
worm surface in is given as : 


= 0 

A^ 1] D^F = 0 


(35) 


the gear surface in S n is given as : 

AjjP&l 2X) F = 0 

Ai^Ai 21) D t F = 0 ( 36 ) 

A e ( . n2) J4 e (21) I3b (1) F = 0 


and the gear surface in is given as : 

A (JC> a£? 2) Aq ZX) F = 0 

A^ n) A^ 21) D^F = 0 ( 37 ) 

A (ja) a£? 2) a£ 21) F = 0 


3.6.2 The Relative Velocity at a Contact Point 

In Eqs . 35 to 37 if one sets 9 = 9 * , the contact lines expressed in the 
stationary systems are obtained for the given 6 position. Then the relative 
velocity at a contact point can be calculated, which is the difference of the 
worm velocity and the gear velocity at this contact point. In the relative 
velocity is given as: 
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( 38 ) 


Vz~Vt-Vt 

v^=-yj (d x I + ^ j 
v£=-y„co 2<-^ ) + -^^>2 > 

where , y ^ , z^ are the coordinates of the contact point in ; Xjj, y n , z n are 
the corresponding coordinates in S n ; i, j and k are the unit vectors in . 

To facilitate further development, let us introduce the term "center 
point" . A center point is the point of intersection between a tooth surface 
and the gear pitch circle. In the meshing zone there is a gear center point 
on every gear tooth. Correspondingly there is a worm center point on the worm 
coil in meshing with this gear tooth. The worm center point is in general 
slightly different than the gear center point, unless the center point is a 
contact point. In the analysis below, the normal and friction forces between 
a pair of meshing teeth are taken as acting at the worm center point, which 
for brevity will be referred to simply as the center point. 

In case that the center point is a contact point, the relative velocity 
is : 


V r = yJvZ+vZ = v /(x i o 1 ) 2 +(d 2 co 2 /2p 
Yj = 0, xj + yl = < <* 2 / 2 > 2 


(39) 


3.6.3 Normal Direction of the Worm Surface 

Because the worm surface is formed by the original generating plane , the 
normal direction at a point of the worm surface is the same as the normal 
direction of the original generating plane. Further, since the worm surface 
and the original generating plane contact in a straight line, all the points 
on the contact line have the same normal direction. In S 2 this normal 
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direction is: 


n x - coSTlr 1 simJr 2 +tanPsini|r 1 
n^-sin^s in\|r 2 +tanp cosilTi 
n z =-cosi|f 2 


(40) 


Transforming Eqs. 40 into S 2 , the normal direction of the gear surface 
at the contact point can also be obtained. Note that since the normal 
direction of the worm is outward from the worm body, the normal direction of 
the gear is then inward to the gear body. The normal direction can be further 
transformed into a stationary coordinate system. Thus, for the given position 
9 * , the normal direction in Sj is: 

n JX = cos (ifri- 6 ^) sini|r 2 + tanpsin 

X 2 j y =-sin(^r x - 8 i) sini|r 2 + tanpcos (41) 

n Jz =-cost|r 2 


3.6.4 Clearance Function 

Once the worm and gear surfaces are found, the clearance between the two 
mating surfaces can be determined. This information is of vital importance in 
the analysis of fluid film lubrication. This clearance is defined to be the 
gap between the two surfaces measured in the direction of the common normal on 
the contact line. Clearly, as the distance from the contact line increases, 
the gap so measured becomes less representative of the true distance between 
the two surfaces. But this is immaterial, because as the gap becomes large, 
the effects of fluid film lubrication also become negligible. The clearance 
function can be obtained numerically as follows: (1) Choose the grid system 

on the gear tooth and obtain the coordinates (x 2 ,y 2 ,z 2 ) °f t * ie points. 

(2) Calculate the normal direction at a contact point and obtain the equation 
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of the straight line passing through a grid point and in this normal 
direction. (3) Find the point of intersection between this straight line and 
the worm surface. The gap at the grid point is then the distance between the 
grid point and the point of intersection. It is convenient to carry out this 
calculation in the lubrication analysis, where the grid system has to be 
selected first. 

3.7 Results and Discussion 

Sample results were computed for the selected wormgear . Fig. 3.6 shows 
the wormgear position at 6 =0, where ten pairs of teeth are in meshing, 
beginning from gear tooth No. 2 and ending with gear tooth No. 11. The contact 
lines on several of these teeth are shown in Figs. 3.7 as viewed in the (x rij 
z n ) plane. In each of these figures the circles are on the gear tooth, whose 
domain is shown by the solid lines; the stars are on that portion of the worm 
thread that overlapps the gear tooth. Note also that the unit of the 
coordinates in these figures is millimeter, and the teeth are shown in their 
true proportion. On gear tooth No. 2, Fig. 3.7(a), the two contact lines are 
far from the mid-plane (z n =0) . These contact lines move toward the mid-plane 
as the tooth number increases. See Figs. 3.7(a) to 3.7(d). On gear tooth 
No. 10, Fig. 3.7(e), the second enveloping contact line crosses over the first 
enveloping contact line. The crossed- over part is nonexistent, as is 
explained previously. On gear tooth No. 11, Fig. 3.7(f), the entire second 
enveloping contact line falls below the first enveloping contact line, and 
becomes non-existent. Thus, toward the end of the meshing zone, the total 
length of the contact lines diminishes . From these figures it can also be 
seen that there are two contact lines on the majority of the teeth in the 
meshing zone; these contact lines are located rather symmetrically on the two 
sides of the mid-plane and they run across the tooth face more or less in the 
direction of the x 2 axis . 
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The total length of the contact lines on each meshing tooth is plotted 
in Fig. 3.8. The situation corresponding to the position 0*=O is shown in 
Fig. 3.8(a). When the gear rotates one degree, i.e. 0*=1 deg, tooth No. 11 
moves out of, while tooth No.l has not moved into, the meshing zone. Hence 
there are only nine pairs of teeth in meshing. The distribution of the 
contact line lengths in this situation is shown in Fig. 3.8(b). When the gear 
rotates one more degree, i.e. 6*—‘l deg, tooth No.l has moved into, while tooth 
No. 10 still remains in, the meshing zone. There are again ten pairs of teeth 
in meshing. The corresponding distribution of the contact line lengths is 
shown in Fig. 3.8(c). Since the angle between two neighboring teeth is 360/N 2 
=3.273 deg, if the gear further rotates 1.27 deg, the situation of Fig. 3.8(a) 
re -appears . 

The variation of the contact line length on one gear tooth as it moves 
through the meshing zone is shown in Fig. 3.9, which is essentially the three 
graphs of Figs. 3.8 combined. 

The radial position of the edge curve on the worm surface is shown in 
Fig. 3.10. For comparison, the radial position of the hour- glass -shaped worm 
pitch cylinder is plotted as the dashed curve. Since the tooth depth is about 
20 mm, it can be seen that the edge curve is located far below the worm 
thread. Hence there is no problem of undercutting the worm thread. 

From the above results one can identify the following merits of the 
selected wormgear: 

1. On the majority of the meshing teeth there are two contact lines. This 
means that the contact region on a tooth of the double -enveloping 
wormgear is larger than that of the single -enveloping wormgear, which 
only possesses one contact line on a meshing tooth. 

2. The contact lines run across the tooth face more or less in the 
direction of the x 2 axis. Since the relative velocity of the meshing 
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surfaces is nearly perpendicular to the x 2 axis, the two directions meet 
at a large angle. This feature is favorable to the establishment of 
fluid film lubrication at the contact. 

3. On the majority of the meshing teeth, the pattern of the contact lines 
is nearly symmetrical with respect to the center point. This feature 
allows a major simplification to be made in the load sharing analysis. 

4. With the chosen gear dimensions there is no problem of undercutting the 
worm thread. 

4 Load Sharing Analysis 

The methodology of determining how the load is distributed among meshing 
teeth is not well developed. Rational analyses are rare and complicated. For 
spur and helical gears, the load sharing analyses are usually based on the 
variation of tooth stiffness along contact lines (Elkholy, 1989). In the 
practical design of these gears, the approach is to analyze an idealized 
meshing model, then apply various factors to modify the results. To estimate 
the magnitudes of these load distribution factors , the common approach is to 
assume that gear teeth behave like a series of thin slices that deflect 
independently (Steward, 1989) . It is doubtful that these methods are suitable 
to the double -enveloping wormgear, because its contact pattern is vastly 
different than those of the spur and helical gears. Hence, for the double- 
enveloping wormgear, there does not exist a theory of load sharing. Neither 
are there adequate experimental data to guide the analysis and design. In the 
following, we shall base our analysis of load sharing on a hypothesis. 

4.1 Hypothesis on Load Sharing 

At a contact point the forces between the meshing surfaces consist of a 
normal component, that acts in the direction of the common normal, and a 
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friction component. The latter acts in the direction of the relative velocity 
and has a magnitude equal to the normal component multiplied by the Coulomb 
friction coefficient. Since the direction of the normal and that of the 
relative velocity vary with position, these forces vary from contact point to 
contact point not only in magnitude but also in direction. However, the 
contact analysis in Section 3 reveals that the contact line pattern on each 
tooth is nearly symmetrical with respect to the center point. In addition, 
because of the large speed ratio (and consequently the small lead angle) of 
the selected wormgear, the variation of the normal direction on a tooth 
surface is small. These two features allow the following simplifications to 
be made : 

1. The resultant normal force on a tooth is proportional to the contact 
line length on the tooth. It acts at the center point and is in the 
normal direction of the center point. 

2 . The resultant friction force on a tooth acts at the center point and is 
in the direction of the relative velocity at the center point. 

The above simplifications may be considered as a hypothesis of load 
sharing, whose validity should be verified later when experimental data or 
results of a deformation analysis become available. 

4.2 Normal Forces on Gear Teeth 

The coordinates (x^.z^) in S i and the parameter if> mv of the worm center 
point for the mth pair of meshing teeth can be found by solving Eqs . 35 and 
the following set of equations describing the gear pitch circle: 

ia-Xj) 2 + z) = ( d z /2 > 2 (42) 

Yj~ 0 

The normal direction at the center point can be calculated from Eqs . 41 , where 
V> mp is to be used for and 6* is the given position. 
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From Fig. 4.1, where friction is not shown, the torque contributed by 
the normal force on the gear tooth of the mth meshing pair is: 


= Fc_ x t 


z e a ~ c ( a- Xp )^i ] 


n. 


n w = Jn 2 jx + n§y+h 


2 

jz 


(43) 


where r em is an effective radius; n^*, nj y , n jz are as given in Eqs. 41; [ ] m 
denotes the mth meshing pair. 


By hypothesis, F ^ is proportional to the length of contact lines. 
Besides, T^STg*.. Hence the formula: 


“ Gm 


L m 

k-1 


(44) 


where 

F^ 11 is the normal force acting on the gear tooth of the mth meshing 
pair, 

T g is the required output torque from the gear shaft, 

T,^ is the length of contact lines on the mth pair of the meshing 
teeth, and 

N p ' is the actual number of meshing pairs. 

The relative motion between the worm and gear is pure sliding. Hence, 
friction plays an important role in the performance of the wormgear (Shigley 
and Mitchell, 1983). Let fi be the Coulomb friction coefficient, then the 
magnitude of friction is /zF gj 1 , and its direction on the gear is the same as 
the direction of the relative velocity v 1 -v 2 . 

At the center point of the mth pair the worm velocity is parallel to the 
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(Xj,yj) plane and equal to x^c^, the gear velocity is in the (Xj,Zj) plane and 
equal to d 2 a> 2 /2. In the direction of the relative velocity is given by 
Eqs . 39 and is shown in Fig. 4.2. 


Let F fGtn n be the normal force on the gear tooth of the mth meshing pair 
when friction is present. The torgue due to the friction force, AiF fGm n , is: 


■ efm 


= [ 


Tfm P'qb z efm 

n 




( 45 ) 


j{x p <a r ) 2 + (d 2 (*> 2 /2) 


2 2 


Then, F fGm n is found to be: 


v n — 
“ fGm 


Tn L„ 


E L ic ( 


( 46 ) 


z ek r efk 


Jc-l 


4 . 3 Tangential Forces on the Worm 


The force acting on the worm can be decomposed into the radial , 
tangential, and axial components, viz. F^ 3 ^ , F^^ , and F^^ (Shigley and 
Mitchell, 1983). See Fig. 4.1. The tangential component balances the input 
worm torque, which in turn determines the input power. The ratio of the input 
power without friction to that with friction, averaged over an interval where 
the gear rotates 3.273 deg, determines the efficiency of the wormgear. 


The tangential force on the worm for the mth pair of the meshing teeth 

is : 

without friction 





( 47 ) 
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with friction 


FfL= Ffcm l 


V X P <*1 

V < *p»x) 2 + ( <* 2<*2 /2) 2 





( 48 ) 


Since the worm torque equals the sum of the torques contributed from the 
meshing teeth, it is: 
without friction 


K 

Fff — yi Fffyj Xjyp 
m-1 


( 49 ) 


with friction 


fw 


~ 5^ FfWm x ap 
m»l 


( 50 ) 


4.4 Results and Discussion 

Load sharing was computed for a power output of 4000 kW and a rotation 
speed of 130 r/min of the gear shaft. To account for the effect of friction 
at the worm and gear contact, a Coulomb friction coefficient value of 0.03 was 
used in the computation. The results are shown in Figs. 4.3 to 4.5. In these 
figures , the stars represent the case with friction and the circles the case 
without friction. 

Figs. 4.3(a) to 4.3(c) present the normal forces acting on the meshing 
gear teeth for the three wormgear positions, viz. 6* = 0, 1 and 2 deg. These 
figures correspond to the cases of Figs. 3.8(a) to 3.8(c). The variation of 
the normal force on one gear tooth as it moves through the meshing zone is 
shown in Fig. 4.4, which is essentially the three graphs of Figs. 4.3 
combined. It is seen that the difference between including and not including 
friction cannot be detected. This implies that, for a given gear output 
power, the lubricant film at the worm and gear contact must withstand almost 
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the same normal force whether friction is present or not. 

The tangential forces acting on the worm thread corresponding to the 
three wormgear positions of 6 — 0, 1 and 2 deg are shown in Figs. 4.5(a) to 

4.5(c). Here, including or not including friction makes a big difference. 
This is because for the selected wormgear the direction of the relative 
velocity, which is the direction of friction, is very close to that of the 
tangential force on the worm. Thus, a large portion of the input torque must 
be spent to overcome the frictional torque. Hence, effective lubrication at 
the contact can drastically reduce the input power required of the worm. 


5 Description of Computer Programs 

Two computer programs were developed in this work, the graphics program 
and the analyses program. The graphics program is for creating the 3-D 
graphics of the selected wormgear. The analyses program is for implementing 
the geometric , kinematic and load sharing analyses . 

5 . 1 Graphics Program 

The graphics program is named NASA.PRG and is written in I-DEAS 
programming language, IDEAL. 

5.1.1 I-DEAS Software 

I-DEAS (Integrated Design Engineering Analysis System) is a widely used 
commercial software (SDRC, 1991). It has three program families: solid 
modeling, finite element and advanced technique. Solid modeling is used to 
geometrically create an object. The tasks in the solid modeling family 
include construction geometry, object modeling, mechanism design, drawing 
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layout , etc . 


Construction geometry is used to create profiles, wireframes and skin 
groups. A profile is defined by a boundary that lies in a plane. The 
boundary can consist of lines, arcs, circles and splines, which can be defined 
by points, centers and radii. In the boundary holes can be defined. 

Object modeling is used to design mechanical parts which are modeled as 
objects. There are several ways to create an object: (1) Use primitives such 

as blocks, cylinders, cones and spheres. (2) Extrude or revolve a profile 
into a solid. (3) Create a group of profiles and then generate a solid 
through the group. This task also provides jointing, cutting and assembling 
options. Once two objects are created, a third object can be created from 
them by using these options. The created object can be displayed with options 
such as scale, hidden line, shade display, etc. The created object can also 
be rotated, moved and oriented to the required position; and saved and 
printed. 

I-DEAS can be run interactively by selecting from the menu on the screen 
using the mouse, or by inputing from the keyboard or through a program. 

I-DEAS has a high level programming capability using the language IDEAL. 
The commands can be written in a program file which is then executed. The 
programming capability permits arithmetic expressions, variables, mathematic 
functions, and algebraic operations; it also provides commands for branching 
and do -looping, as well as controlling the input and output. 

5.1.2 Description of the Graphics Program 

The process of creating a double -enveloping wormgear by using I-DEAS is 
similar to that of manufacturing the wormgear. This is carried out as 
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follows : 


First a worm blank and a disc grinding wheel are created. Once these 
two objects are created, the worm can be generated by cutting the worm blank 
with the grinding wheel. In the worm creating process the motion of the worm 
blank and the grinding wheel is similar to that of the first enveloping 
process . 

Then the gear blank is created. The gear is created by cutting the gear 
blank with the previously created worm. The motion of the gear blank and the 
worm is similar to that of the second enveloping process. 

The motions are simulated by rotating the blanks and cutters about their 
respective axes. The rotation is divided into incremental steps. At each 
step the cutting operation takes place. The incremental step determines how 
smooth the generated worm and gear surfaces are. As a result there is a 
difference between the I-DEAS generated surface and the manufactured surface. 
Because of the computer storage limitation, the steps can not be chosen too 
small . 


The structure of the graphics program is shown by the flow chart in Fig. 

5.1, which consists of the following steps: 

1. Input wormgear dimensions and parameters which are needed in its 
manufacture . 

2. Input the incremental stepsizes ('dtw' for generating the worm and 'dtg' 
for generating the gear) and the total number of steps in the cutting 
processes ('wn' for generating the worm and 'gn' for generating the 
gear) . 

3. Construct the grinding wheels for generating the worm. Four are 
constructed and jointed in sequence to reduce the CPU time, i.e. four 
coils of the worm thread are created for one revolution of the worm 
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blank (dtw*wn=360 deg) . 

4. Create the worm and gear blanks. 

5. Create the worm. When the process is completed, the created worm is 
saved under the name 'worm' . 

6. Create the gear. When the process is completed, the created gear is 
saved under the name 'gear'. 

7. Display the created worm and gear. 

The program is listed in Appendix 1. The program works interactively 
and ends at the object modeling screen. This allows the user to see the 
created objects again before exiting I-DEAS. 

The procedure to execute the program is as follows: (1) Enter I-DEAS and 
come to the construction geometry task screen. (2) Run program in the 
construction geometry screen. More detailed instructions are given in 
Appendix 1 . 

Because the worm thread of the double -enveloping wormgear is asymmetric 
about its center, the number of steps needed to create the whole worm is 
larger than that for creating other machine parts , such as the spur gear or 
the cylindrical worm. More coils of the worm thread requires more steps. 

Every step increases the number of surfaces on the worm blank. The larger the 
number of surfaces , the more CPU time and computer disc space are needed to 
calculate and save the data. Once the capacity of the computer or the I-DEAS 
space is full, the screen stops, or the program crashes. This is why there 
are only four coils of worm thread shown in Fig. 2.7. The program was 
executed on an IBM 5080 terminal at SUNY- Binghamton. 

From the above description it can be seen that (1) the computer 
simulation is helpful to visualize the selected wormgear, and (2) the computer 
simulation can not replace the analysis, because the numerical data of the 
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generated surfaces are not accurate enough for the projected lubrication 
analysis . 


5 . 2 Analyses Program 


5.2.1 Overall view of the Program 


The analyses program is named NASA. FOR and is written in FORTRAN. This 
program consists of a main program and several subroutines. The main program 
is the driver program of all the subroutines and may be tailored and adapted 
to the final design package. The subroutines fall in two categories: those 
driven by the main program (DIMEN, WORM, WORML, GEAR, CONTACT, PITCH and LOAD) 
and those driven by other subroutines (EQW, EQWL, EQG, FIRST, SECOND, WORMT , 
ROOT and NORMAL). Most of the logical judgements are done in the subroutines 
called by the main program. The evaluations of analytical expressions are 
done in the subroutines called by other subroutines. 

The structure of the analyses program is shown by the flow chart in Fig. 
5.2, which consists of the following steps: 

1. Calculate the wormgear dimensions (DIMEN). 

2. Calculate the coordinates in S x of the worm surface and the curvatures 
of the worm surface (WORM and WORML) . 

3. Calculate the coordinates in S 2 of the points on the gear surface along 
the contact line and the curvatures of the gear surface (GEAR) . 

4. Calculate the contact lines in S n on every pair of the meshing teeth for 
a given position (CONTACT) . 

5 . Calculate the coordinates in and the kinematic properties of the worm 
center points for a given position (PITCH) . 

6. Calculate the load distribution on the meshing teeth for a given 
position (LOAD) . 
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Input of this program, under the ' parameter ' command, are at the 
beginning of the main program and are explained there. In the main program 
the given position 8 * is named 'tn2', which should be the input for the 
subroutines CONTACT, PITCH and LOAD. 

Input to a subroutine are explained at the beginning of the subroutine. 
The output to be used by other subroutines are described in the subroutine 
statements . The output not to be used by other subroutines are stored in data 
files (DIMEN, WORM, WORML, GEAR, CONTC, PITCH and LOAD). 

Error flags would be shown on the screen such as 'dimension is not big 
enough', 'the stepsize is too large', 'the worm surface intersects the edge 
curve' , etc . 

5.2.2 Detailed Description of the Program 
Main Program 

The input under the 'parameters' command include the center distances a 0 
and a, the speed ratios i 0 and i, the number of worm thread N x , the number of 
meshing teeth N p , the worm pitch diameter d ls the apex angle £, the rotation 
speed of the worm, the output power of the wormgear, and the friction 
coefficient at the contact. 

The stepsizes and the dimension arrays are input under the 'data' 
command. The input stepsizes are: 'df' (for the running parameter Vh) > '<*t' 
(for the running parameter 8 2 ) , ' dx' (in the x axis), and 'dz' (in the z 
axis). There are four one -dimensional common arrays with dimension nn=600. 

The given position 8 * is named 'tn2' . The wormgear position of 8 =0 is 
shown in Fig. 3.6. 'tn2' is not to be larger than the angle between two 
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neighboring gear teeth. The two angles 't20' and 't21' specify the interval 
of the running parameter 0 2 , which is chosen to be slightly larger than the 
meshing zone . 

The convergence criteria for root search are specified by 'ep' (for 
angles) and 'epp' (for distances). A logical judgement for the number of 
meshing teeth is set up in the main program. This is because the actual 
number of meshing teeth obtained from subroutines CONTACT and PITCH should 
agree, otherwise the results can not be used in the subroutine LOAD. A flow 
chart for the main program is shown in Fig. 5.3. 

Subroutine DIMEN 


This subroutine must be run before any other subroutine, because it 
provides the wormgear parameters and dimensions. The formulas used for the 
calculation can be found in AGMA (1965) , Dudley (1984) , and Shigley and 
Mitchell (1983) . 

The meshing zone is specified by 'start' and 'end', which are the angles 
where a gear tooth runs into and out of the meshing zone respectively. 'f20' 
and 'f21' are the two angles specifying the interval of the running parameter 
V> 2 . 'tao' is the angle between two neighboring gear teeth. The half gear 
face angle, named 'dlt' , is used to specify the Vi interval in the calculation 
of the gear surface. 'reps' is an angle for determining the position of the 
original generating plane that grinds the other side of the worm thread, Fig. 

5.4. 

Subroutine WORM and WORML 

These subroutines implement the analysis described in Sections 3.3 and 

3.5. The two sides of the worm thread are calculated by WORM and WORML, 
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respectively. 


Three logical judgements are set up. The first checks if the worm 
surface intersects the edge curve. The second checks which branch of the 
envelope surface is the real worm surface, Fig. 3.2. The third examines the 
tooth domain. 

Subroutines EQW and EQWL are called by subroutines WORM and WORML, 
respectively. A flow chart for WORM (and WORML) is shown in Fig. 5.5. 

Subroutine GEAR 

This subroutine implements the analysis described in Sections 3.4 and 
3.5. Calculations for both the conventional form and the modified form are 
provided. 

A 6 value, named 'tl', is first given. Then ip y named 'fl', is varied 
near the 6 value to calculate the coordinates of the points on the gear 
surface. This is because if ip is far away from 6, it is also off the tooth 
domain. The range of ip variation is twice the gear face angle. This range is 
much less than the range of the running parameter ip in the subroutine WORM. 
Hence, the stepsize Aip lt 'df', used here is a tenth of the stepsize given in 
the main program. 

The logical judgements set up in this subroutine consist of the tooth 
domain examinations and sorting out the fictitious points. A flow chart for 
GEAR is shown in Fig. 5.6. 

Subroutine CONTACT 


This subroutine implements the analysis of Section 3.6.1. Calculations 
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for both the conventional form and the modified form are provided. At a given 
position 'tn2', the contact lines on all the meshing teeth are calculated. 

The angles 'start' and 'end' are used in this subroutine to find out how many 
pairs of meshing teeth actually exist in the meshing zone. 

Similar to the case of GEAR, ip is varied near the 6 value. To obtain 
accurately the length of the contact lines, as many points are needed as is 
practicable. Hence the stepsizes of if) and z are chosen to be df/100 and 
dz/10 , where 'df' and 'dz' are the input stepsizes from the main program. 

Four temporary one - dimens ional arrays are used to save the coordinates 
of the points on the same contact line. 'nn' gives the dimension of these 
arrays, which is controlled by the main program. If the saved points exceed 
'nn' , an error flag would appear on the computer screen. The number 'nc' 
(output name 'ncc') records how many pairs of meshing teeth exist in the 
meshing zone for the given position, which should be the same as that obtained 
from the subroutine PITCH. 

Output of the contact line length is stored in the array 'CONTC', which 
has dimension N p and is the input of the subroutine LOAD. A flow chart for 
CONTACT is shown in Fig. 5.7. 

Subroutine PITCH 


This subroutine implements the analysis of Sections 3.4 and 3.6.2. At 
the given position 'tn2' (named 'fn2' here), the worm center point on each 
pair of the meshing teeth are obtained by first finding the points on the worm 
profile in the mid-plane, i.e. y^ = 0, then sorting out the center point on 
the profile by using the subroutine ROOT. The angles 'start' and 'end' are 
used to check if a worm coil falls in the meshing zone. 
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To obtain accurately the center point coordinates while still saving 
the CPU time, the stepsize is chosen to be df/100, where 'df' is the input 
steps ize from the main program. The range of if> variation is 90 deg around the 
given position 'fn2'. The index mw=l would mean that the worm profile 
intersects the edge curve. The coordinates of the center points obtained are 
in Sj. These coordinates and the normal directions at the center points are 
stored in one -dimensional arrays with the maximum dimension N p , and are input 
for the subroutine LOAD. 

The number 'nc' (output name 'ncp') records how many pairs of meshing 
teeth exist in the meshing zone for the given position, which should be the 
same as that obtained from the subroutine CONTACT. A flow chart for PITCH is 
shown in Fig. 5.8. 

Subroutine LOAD 


This subroutine implements the analysis of Sections 4. Before this 
subroutine is run, the subroutines PITCH and CONTACT must be run so as to 
provide the input for this subroutine. 

Subroutine ROOT 


This subroutine searchs the root in a given interval where only one root 
exists. The root searching is based on the idea of bisection (Press, 1986), 
but three sections are used instead to reduce the CPU time . Here the root is 
a point at which the radius of the worm is equal to the radius of the worm 
pitch cylinder at the same axial position. The criteria used to search the 
center point are the difference of the two radii or the angles between the two 
grid points. 
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6 Summary 

Based on the consideration of achieving better contact properties 
(longer contact lines and a larger angle between the contact line and the 
relative sliding velocity), a double -enveloping wormgear was selected. The 
worm was generated by a disk grinding wheel. The center distance and the 
speed ratio were respectively the same in the first and the second enveloping 
processes (conventional form) . 

Based on the suggested speed ratio and center distance for helicopter 
applications , the other dimensions and geometric parameters of the wormgear 
were laid out basically in accordance with the AGMA gear standards. However, 
some of these quantities were modified as a result of compromising between 
conflicting requirements . It is felt that for a novel application like ours 
that is not covered by past experiences, some deviations from the common 
practice would be inevitable. 

3-D graphics representation of the selected wormgear was developed using 
the I-DEAS software. The software works like a cutting tool generating a 
surface from a blank, but in discrete steps. Because of the computer storage 
limitation, the steps could not be made too small. As a result, the numerical 
data of the generated surfaces were not accurate enough for the projected 
lubrication analysis. However, the graphics representation does convey a 
clear notion of the proportion of the selected wormgear. 

The kinematics of meshing was analyzed with the method of enveloping. 

The method is simple and elegant, and can be used to construct conjugate 
surfaces in other contacting mechanisms as well. Using this method the shape 
of the worm and gear surfaces , the contact lines , the normal directions , the 
relative normal curvatures , and the relative sliding velocities at the contact 
points were obtained. With these results the clearance function between the 
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meshing surfaces can also be obtained, which is needed for the projected 
lubrication analysis. 

Based on the contact line pattern revealed by the kinematics analysis, a 
hypothesis of load sharing was made. Using this hypothesis the variation of 
the normal contact force acting on a gear tooth as it ran through the meshing 
zone, as well as the distribution of the normal forces on all the meshing gear 
teeth at any instant, were obtained. This information is needed in the 
projected lubrication analysis. The required torque output of the worm was 
likewise obtained. The results clearly demonstrated the advantage and need of 
implementing effective lubrication at the contact. 

A computer program, written in FORTRAN, was developed to carry out the 
computations. The structure of the computer program and its usage are 
described in detail. 
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Table 1 Dimensions of the Selected Wormgear 


TERM SYMBOL 

VALUE 

FORMULA BASED 

REFERENCE 

speed ratio 

i 

no 



given 

normal pressure 
angle 

Q n 

20 deg 



Dudley, 1984 
p. 3.71 

center distance 

a 

880 mm 



selected 

number of worm 
threads 

Ni 

1 



selected 

number of gear 
teeth 

n 2 

no 

n 2 = 

N x *i 


pitch diameter 
of worm* 

di 

300 mm 

d, - a 0 875 Aa (E) 

(k d - 1.7-2. 2) 

Dudley, 1984 
p. 3.70 

pitch diameter 
of gear 

d 2 

1460 mm 

d 2 — 

2*a-d x 


axial pitch 

Px 

41.6975 mm 

Px = 

7T*d 2 /N 2 

Dudley, 1984 
Table 3.32 

worm lead angle* 

7 

2.533 deg 

tan7 

= Px* N iA/di 

Dudley, 1984 
Eq. 3.38 

normal cirlular 
pitch 

Pn 

41.6568 mm 

Pn - 

p x *cos 7 

Dudley, 1984 
Table 3.32 

axial pressure 
angle* 

<*x 

20.018 deg 

tana. 

K = tana n /cos 7 

Dudley, 1984 
Eq. 3.40 


Aa 

0.818 deg 

Aa = 

Px/2/d z 

Dudley, 1984 
p. 3.72 

base circle 
diameter 

d b 

519.3182 mm 

d b = 

d 2 *sin(a x +Aa) 

Dudley, 1984 
Eq. 3.41 

module 

Xm 

13.2727 
mm/ tooth 


Px/^ 

Dudley, 1984 
Table 3.32 

number of 
meshing teeth 

N P 

10 



selected 

half angle of 
meshing 


15.627 deg 

= 

7T*(N p -0.45)/N 2 

Shen, 1983 
Table 10.1 

start angle 
of meshing 

i> f 

5.209 deg 

V» f - 

sin _1 (d b /d 2 )--0 a 

Shen, 1983 
Table 10.1 

thickness of 
worm thread 

Si 

18.7639 mm 

Si = 

0 . 45*p x 

Dudley, 1984 
p. 3.73 

thickness of 
gear tooth 

s 2 

22.9336 mm 

s 2 = 

0 . 55*p x 

Dudley, 1984 
P. 3.73 
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whole depth 

b t 

20.8284 mm 

working depth 

b* 

18.7455 mm 

addendum 

b a 

9.3728 mm 

clearance 

c 

2.0829 mm 

throat diameter 
of worm* 

dal 

318.7455 mm 

root diameter 
of worm* 

dfi 

277.0888 mm 

throat diameter 
of gear 

d a 2 

1478.7455 mm 

root diameter 
of gear 

d f2 

1437.0888 mm 

face width of 
worm 

b i 

393 mm 

face width of 
gear 

b 2 

249 mm 

apex angle of 
generating plane 


5 deg 


b t - 0.5*p n 

Dudley, 1984 
Table 3.31 

b k = 0.45*p n 

Dudley, 1984 
Table 3.31 

b a - 0 . 225*p n 

Dudley, 1984 
Table 3.31 

c - b t -b k 

d a i - d x +2*b a 

Shen, 1983 
Table 10.1 

dfi = d al -2*b t 

Dudley, 1984 
Eq. 3.39 

d a 2 “ d 2 +2*b a 

Dudley, 1984 
Table 3.32 

df2 = d a2 -2*b t 

Shen, 1983 
Table 10.1 

b x = d 2 *sinV' a 

Shen, 1983 
Table 10.1 

b 2 = (0.9-1. 0)*d fl 

Shen, 1983 
Table 10 . 1 

selected 


* At the center of worm. 
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Figures 


Fig. 2.1 Single -enveloping wormgear (cylindrical wormgear) 

Fig. 2.2 Double -enveloping wormgear 

Fig. 2.3 Worm ground by a disc grinding wheel 

Fig. 2.4 Worm ground by a cone grinding wheel 

Fig. 2.5 Double -enveloping wormgear layout and nomenclature 

Fig. 2.6 Selected wormgear 

Fig. 2.7 Selected worm 

Fig. 2.8 Selected gear 

Fig. 2.9 Close-up view of the selected gear teeth 

Fig. 3.1 Coordinate systems 

Fig. 3.2 Envelope surface (worm) and its edge curve 

Fig. 3.3 Contact lines on a gear tooth of the modified form 

Fig. 3.4 Contact lines on a gear tooth of the conventional form 

Fig. 3.5 Worm and gear tooth domains 

Fig. 3.6 Wormgear position at 6* = 0 

Fig. 3.7 Contact lines on the meshing gear teeth 

Fig. 3.8 Length of contact lines on the meshing teeth 

(a) position of 0*-O deg 

(b) position of 8*-l deg 

(c) position of 6 =2 deg 

Fig. 3.9 Variation of contact line length on a gear tooth 
Fig. 3.10 Radial position of the edge curve 
Fig. 4.1 Normal force and its components 
Fig. 4.2 Friction force and its components 

Fig. 4.3 Normal forces on the gear teeth 

(a) position of 0*=O deg 

(b) position of 0*=1 deg 

(c) position of 2 deg 

Fig. 4.4 Variation of normal force on a gear tooth 

Fig. 4.5 Tangential forces on the worm 

(a) position of 0*=O deg 

(b) position of 0*=1 deg 

(c) position of 6 =2 deg 

Fig. 5.1 Structure of graphics program 
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Fig. 5. 
Fig. 5. 
Fig. 5. 
Fig. 5. 
Fig. 5. 
Fig. 5. 
Fig. 5. 


2 Structure of analyses program 

3 Flow chart of main program 

4 Two sides of the worm thread 

5 Flow chart of subroutine WORM (also WORML) 

6 Flow chart of subroutine GEAR 

7 Flow chart of subroutine CONTACT 

8 Flow chart of subroutine PITCH 
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Fig. 2.3 Worm ground by a disc grinding wheel 
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Fig. 2.5 Double-enveloping wormgear layout and nomenclature 
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Fig. 2.6 Selected wormgear 
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Fig. 2.8 Selected gear 
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Fig. 2.9 Close-up view of the selected gear teeth 
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Fig. 3.2 Envelope surface (worm) and its edge curve 
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Fig. 3.3 Contact lines on a gear tooth of the modified form 
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Fig. 3.7 Contact lines on the meshing gear teeth (c) 
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Fig. 3.8 Length of contact lines on the meshing teeth 
(a) position of 6*=0 deg 










©Ajno 0Bp9 Qi\\ p uojijsod |e|pey (H'0 ‘6jd 


(««n) -Z 



Radius (mm) 



Fig. 4.1 Norma! force and its components 
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Fig. 4.2 Friction force and its components 
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Fig. 4.3 Normal forces on the gear teeth 
(a) position of 0*=O deg 
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Fig. 4.3 Normal forces on the gear teeth 
(b) position of B*=1 deg 
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Fig. 4.4 Variation of normal force on a gear tooth 
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Fig. 4.5 Tangential forces^ on the worm 
(b) position of 0*=1 deg 
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Fig. 5.1 Structure of graphics program 


80 










Fig. 5.2 Structure of analyses program 
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Call LOAD 

calculate the normal load 
and the tangential component. 



Fig. 5.3 Flow chart of main program 
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Fig. 5.7 Flow chart of subroutine CONTACT 
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Fig. 5.7 (continued) 
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PITCH 


Input 


ind=0, nc=0 


m=(t21-t20)/dt 

n=(pi/2.0)/df 

yj”° 


j=l,n 


Call WORMT 


-""Ts the 

point on the edge 
^curve ? 


As the point on ^ 
'-the right branch ? 

^ls the point in the. 
tooth domain 


ind=ind+l 
r(ind)=(x, z) 
f(ind)= fl 

<i>^ 

{Y es 


k=l, ind 

Find interval of iir x 
where the center point on 



Call ROOT 
find the center point 
Oi, x, z) 



Fig. 5.8 Flow chart of subroutine PITCH 
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Call NORMAL 
calculate normal direction 
and relative velocity. 



Fig. 5.8 (continued) 
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Appendix 1 Graphics Program (NASA.PRG) 

The procedure to execute r NASA.PRG f 

(1) From the SUNY - Binghamton operating system prompt (assumed here to be $) : 
$ ready: type 'runideas' <return> 

$ jump to IBM screen (use Alt- jump) 

$ user i.d.: type 'ideas2' 

$ password: type 'no' 

$ jump back to BINGVMB screen type 'none' (use Alt-jump) 

$ disk space : type a number such as ' 15000 ' 

(waiting for temporary disk to be created) 

$ user i.d.: type your own i.d. number and read password 
$ terminal type: type '15080' 

$ jump screen (use Alt-jump) 

(waiting for I-DEAS screen to appear) 

(2) Under the I-DEAS screen, from the I-DEAS operating prompt (assumed here 
to be #) : 

Enter model file name (none)# <return> 

New model file? (yes)# <return> 

Enter model description (none)# <return> 

Enter system of units (SI)# <return> 

Enter operation mode (Normal_Operation)# <return> 

(3) Then choose which “family" of I-DEAS you want. Notice that choices are 
in the form of a 'menu' . So you can pick a choice with the mouse or type 
under the I-DEAS operating prompt (assumed here to be #) : 

Select menu# SO <return> 

(waiting for Solid Modeling Screen to appear) 

Select menu# CG <return> 
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Select menu# MF <return> 

Select menu# PR <return> 

Select menu# R <return> 

File name# NASA <return> 

Then the program is run interactively. All the results will be seen from the 
screen. 

(4) The created object can be observed again before exiting I-DEAS system. 
The procedure to observe the created object is: 

Select menu# HA <return> 

Select menu# G <return> 

Enter object name# type object name (worm or gear) 

(the created object appears on the screen) 

Once the created object is obtained, choose from the menu options such as 
orient, view, display and AU (auto scale) to observe the object. 

(5) When the observation is completed, put the created object away from the 
screen by: 

Select menu# / <return> 

Select menu # HA <return> 

Select menu# PU <return> 

(6) When all the observations are completed, if needed the created object 
can be saved by: 

Select menu # / <return> 

Select menu# HA <return> 

Select menu# WR <return> 

Select menu# U <return> 

Enter the file name# type the name you choose <return> 
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<return> 


(7) Exit I-DEAS screen by: 

Select menu# / <return> 
Select menu# type 'exit' 
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Program Listing 


C : 

C : 

C : 

C : 

C : 

C : 

C : 

C : 

C : 

C : 

C : 

C : 

C : 

C : 
C : 
C : 
C : 
C : 
C : 
C : 
C : 
C : 
C : 
C : 
C : 
C : 
C : 
C : 
C : 
C : 
C : 
C : 
C : 
C : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
C : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
C : 


****************************************************************** 
This is an ideas program to construct a wormgear set by the 
principle of enevlop. The constructing process is same as 
the manufacturing process exactly. To cut four worm threads 
the 'wn' must be equal to 360/dtw and 'gn' must be 4*360/dtg. 


The basic parameters provided are : 

a 0 center distance of cutting (a0=880.0 mm) 

db base circle diameter (db— 519.3182 mm) 

bl face width of worm (bl=383.0 mm) 

b2 face width of gear (b2=249.0 mm) 

d2 pitch diameter of gear (d2=1460.00 mm) 

da2 throad diameter of gear (da2=1478 . 7455 mm) 

df2 root diameter of gear (df 2=1437 .0888 mm) 

dfl root diameter of worm (dfl=277 . 0888 mm) 

ratio- -speed ratio of wormgear (ratio=110) 

n2 number of gear teeth (n2=110) 

dtl angle of worm turning in each cut 

dt2 angle of gear turning in each cut 

wn times of worm cutting 

gn times of gear cutting 

bta apex angle of generating plane (bta=10 deg) 

psif angle of meshing start (psif=5.209 deg) 

re p S angle of two cutting surfaces (reps=39 . 8724 deg) 

** 'reps' is calculated by the following formular: 
reps=180-{2*[90-arcsin(db/d2)]+s2/r2) , here 

s2=0.55*pi*d2/n2 i.e. s2/r2=l.l*180/n2 (deg) 
****************************************************************** 


= basic parameters input = 


//echo list 
#a0=880.0 
#n2=110 
# ratio=110 
#db=519 . 3182 
#bta=5 . 00 
#dl=300 . 0 
#d2=1460 . 0 
#dal=318 . 7455 
#da2=1478 . 7455 
#dfl=277 .0888 
#df 2=1437. 0888 
#b 1=393 .0 
#b2=249.0 

#rb=db/2 . 0 
#rl=dl/2 . 0 
#r2=d2/2.0 
#ral=dal/2.0 
#ra2=da2/2.0 
#rfl=df 1/2.0 
#rf 2=df 2/2 . 0 
#dh=2.0829 
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K : #psif=5 . 209 
K : #reps=*39 . 8724 

K : #aa=a0 - sqrt ( ( rf 2+dh) **2 -bl*bl/ 4.0) 

K : #bb=sqrt(ra2*ra2-rb*rb) 

K : #cc=a0-sqrt(rl*rl-b2*b2/4.0) 

K : #dd=a0-sqrt(ra2*ra2-bl*bl/4.0) 

K : #ee= a0-ra2 
K : #dtw=8 . 00 
K : #dtg=15.0 
K : #wn=5 
K : #gn=3 
C : 

c ; — 

C : = creat 2-D profile of worm blank — 

C : =—===——— —=—^^ 

C : 

K : "CR" 

K : "PR" 

K : "K" 

K : 0 ,r2 
K : a0-cc,r2 

jy . 11 F ,f 

K : aO-cc ,bl/2 . 0 

. Itpil 

K : aa,bl/2 . 0 

^ . npn 

K : a0-rf2-dh,0 
K : aa, -bl/2 . 0 

J ^ • It p It 

K : aO - cc , -bl/2 . 0 
K : ”F" 

K : aO-cc, -r2 
K : "F" 

K : 0 , -r2 
K : "F" 

K : "C" 

K : V" 

K : "MA" 

K : "PR" 

K : "STO” 

K : pi 
K : "/" 

K : "DE“ 

K : "P" 

K : "Y" 

K : "/" 

C : 

C : 

C : = creat 2-D profile of worm root blank - 

C : 

C : 

K : "CR" 

K : "PR" 

K * "K" 

K : 0,bl/2 
K : dd,bl/2 
K : "F" 

K : ee , 0 
K : dd, -bl/2 . 0 

JT . flptt 

K : 0 , -bl/2 .0 
K : "F" 
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K 

"SPH" 


K 

"A" 


K 

"C" 


K 

6 


K 

"D" 


K 

ra2+2*dh 


K 

•i 


K 

M CO" 


K 

ir p ii 


K 

n*pii 


K 

"K" 


K 

- rb , 0 , 0 


K 

- rb , aO , 0 
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K : 

K : 

K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
C : 
C : 
C : 
C : 
C : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 


0,0,rb/tan(bta) 


ft pit 
ii*ptr 

"K" 

rb*cos(reps) , -rb*sin(reps) ,0 
rb/cos ( reps ) , 0 , 0 
0,0, rb/tan(bta) 

"N" 


"OR” 

"MO" 

"D" 


0 , 1,0 

aO 

"R" 

"K" 

0 , aO , 0 
0 , 0 ,psif 

it y'ti 

"MA" 

"STO" 

c2 


= creat four connected 3-D worm cutters = 


n ^ I? 

"OR" 

"R" 

11 K M 

0 , aO , 0 
0,0, 360/n2 

H yn 

"MA" 

"G" 

c2 

ii /n 

"CO" 

•ijii 


n yn 

"MA" 

"STO" 

c2 

ii Y« 
n jn 

"OR" 

"Rll 

"K" 

0 , aO , 0 

0,0, 2*360/n2 

H y ii 

"MA" 

" G" 
c2 

ti yn 

"CO" 

"J" 
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K : 

K : 

K : 

K : 

K : 

K : 

C : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 

C : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
C : 
C : 
C : 
C : 
C : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 


ll pH 

f| Kj* 

0,0, -b2/2 
"D" 

0 , 0,1 

ivpti 

lip If 
•I p II 

I* K ,f 

0,0,b2/2 

"D" 

0 , 0,1 

"N" 

"P«* 

II p It 
IIJTII 

0 , aa+100 , 0 
"D" 

0 , 1,0 

"N" 

II ^ II 

"MA" 

,'STO" 

c2 

iiyti 
it yit 

’’MO" 

"DE" 

tty” 

•t j it 

"MA'* 

"PU" 

c2 

ityit 


= creat 3-D worm blank = 



"DEF" 


"R" 

"X" 

32 

"F" 

’’A" 

11 

"D" 

"D" 

"REV" 

Pi 

It yMl 

"OR" 

"R" 

0,0,90 

tty/ it 

"MA" 
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K : "STO" 

K : wl 
K : "G" 

K : c2 
C : 

C : 

C : = 3-D worm cutting = 

C : 

C : 

K : "/" 

K : "CO” 

K : "CU" 

K : 

K : "/" 

K : #time=l 
K : #start : 

K : #tl=time*dtw 
K : #t2=tl/ratio 
K : "OR" 

K : "R" 

K : 

K : - dtw ,0,0 
K : 

K : "HA" 

K : "G" 

K : c2 
K : 

K : "OR" 

K : "R" 

K * ,, K n 
K : 0 , aO , 0 
K : 0,0, t2 
K : "/" 

K : "CO" 

K : "CU" 

K : 

K : "/" 

K : #time=time+l 

K : #if (time It wn) then goto start 
K : "/" 

K : "OR" 

K : "R" 

K : 

K : ( time -l)*dtl, 0,0 
K : 

K : "MA" 

K : "STO" 

K : worm 
K : 

K : "/" 

K : "MO" 

K : "DE" 

K : "Y" 

K : 

K : "CR" 

K : "REV" 

K : p2 
K : 

K : 

K : "OR" 

K : "R" 

K : 
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K : 0,0,90 
K : "/" 

K : "CO" 

K : "J" 

K : 

K : "/" 

K : "MA" 

K : "PU" 

K : worm 
K : "Y" 

C : 

C : ■— 

C : = 3-D gear cutting = 

C : 

C : 

K : "/" 

K : "MA" 

K : "G" 

K : g2 
K : "G” 

K : worm 
K : "/" 

K : "CO” 

K : "CU" 

K : 

K : "/" 

K : #time=l 
K : //again : 

K : #tl=time*dtg 
K : // dt2=dtg/ratio 
K : "OR" 

K : "R" 

K * "K" 

K : 0 , aO , 0 
K : 0,0, dt2 
K : "/" 

K : "MA" 

K : "G" 

K : worm 
K : "/" 

K : "OR" 

K : "R" 

K : 

K : -tl,0,0 
K : "/" 

K : "CO" 

K : "CU" 

K : 

K : "/" 

K : //time=time+l 

K : //if (time It gn) then goto again 
K : "/" 

K : "OR" 

K : "R" 

K * «k m 
K : 0 , aO , 0 

K : 0 , 0 , - (time-l)*dt2 
K : "/" 

K : "MA" 

K : "STO" 

K : gear 
K : 
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K : 7" 

K : "MO” 
K : "DE" 
K : "Y" 

K : "/" 

K : "MA" 
K : "PU" 


K : 
K : 
K : 
C : 
C : 
C : 
C : 
C : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
K : 
C : 
C : 
C : 
C : 
C : 
K : 
K : 
K : 


gear 

MY” 

//echo all 


= view worm = 


"G" 

worm 

ii yn 
tiyii 
ft £ li 

1,1,1,60 

"MG" 

IITMI 

"MG" 

«« t *• 

“DO” 

"HI" 

II £.1 
II ^11 

"DE" 

ii y*» 


= view gear = 


II 

"MA" 

"G" 


K : gear 
K : "/" 

K : "V" 

K : "RM" 

K : 0,0,120 
K : "/" 

K : "DO" 

K : "SH" 

K : "E" 

C : 

E : **** END OF SESSION **** 
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Appendix 2 Listing of the Analyses Program (NASA. FOR) 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


******************************************************************** 
PROGRAM MAIN: 

This is a part of a developing software package which is used to 
analyse the double -enveloping wormgear fluid film lubrication. 

This part includes: 

(1) calculation of the wormgear dimensions; 

(2) calculation of the coordinates of worm and gear surfaces 
in SI and S2 systems; 

(3) calculation of the contact lines on all pairs of meshing 
teeth for a given position; 

(4) calculation of the coordinates and relative velocity at 

the worm pitch points on all pairs of meshing teeth for 
a given position; __ 

(5) calculation of load sharing among all meshing teeth for 
a given position. 


To 


run this 
aO 
a2 
bta 
dl 
dfl 
dt2 
dx 
dz 

ep , epp 

nl 

np 

nn 

rfn 

ratio 

ratio2 

power 

rpm 

tn2 

uf 


program the input parameters are: 
center distance of the first enveloping 
center distance of wormgear 
apex angle of the generating plane 
pitch diameter of worm 

51 running step size in the first enveloping 

52 running step size in the second enveloping 
interval along x axis 

interval along z axis 

convergence criterion for root searching 

number of worm thread 

number of meshing teeth 

number of dimension 

normal presure angle 

speed ratio of the first enveloping 

speed ratio of wormgear 

output power from gear shaft 

input speed from worm shaft 

f iven gear position 
riction coefficient 


(mm) 

(mm) 

(deg) 

(mm) 

(deg) 

(deg) 

(mm) 

(mm) 


(deg) 


(kw) 

(rpm) 

(deg) 


The variables in the main program are: 

t20 running angle of the second enveloping start (rad) 

t21 running angle of the second enveloping end (rad) 


Subroutines included: 


(1) Called by MAIN Program 


DIMEN 

WORM 

WORML 

GEAR 

CONTACT 

PITCH 


develop the dimensions of wormgear. 

develop the worm surface in SI. Here the worm is 

ground by a disc grinding wheel. 

develop the worm surface on the other side, in Si. 
develop the gear surface in S2. 

calculate the length of contact lines on each pair 

of meshing teeth at the given position. 

finding the pitch points on the worm surface, the 
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LOAD 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


c 


c 

c 

c 

c 

c 


c 


normal directions and relative velocities on each 
pair of meshing teeth for a given position, 
calculate the normal force on each pair of meshing 
teeth and the torque sharing among the worm meshing 
teeth. 


(2) Called by Subroutine 

EQG the gear surface equation in S2 

EQW the worm surface equation in SI 

EQWL the worm surface equation of the other side in SI 

FIRST calculation of the first enveloping surface 

NORMAL calculation of the properties of the worm pitch points 
ROOT finding the pitch point 

SECOND calculation of the second enve lop ing^ surface 

WORMT calculation of the worm profile on mid-plane in Sj 


Outputs of the program are written in the data files under 
the names of the subroutines. Units used in the output data are: 

rad for angle 

mm for length 

newton for force. 


PROGRAM MAIN 

implicit double precision(a-h, o-z) 

parameter (bta=5 . 0 , a0=880 . 00 , a2=880 . 00 , ratio=110 . , ratio2=110 , 
& rfn=20 ,nl=l ,np=10 ,dl=300 . 0 , epp=l . 0e-5 , ep=l .0e-3 , 

& power=4000 , rpm=14280 ,uf=0 . 03 ,nn-600) 

dimension contc(np) , tp2(np) , fp2(np) ,pxj (np) ,pzj (np) , 

& xpn(np) ,ypn(np) ,zpn(np) ,vx(np) ,vy(np) ,vz(np) , 

& wtor(np) ,wtorf(np) , forn(np) ,fornf(np) , 

6c exaa(nn) , exbb (nn) , excc (nn) , exdd (nn) 

data df 1/10 . 0/ , dt2/l . 0/ , dx/5 . 0/ , dz/5 . 0/ , tn2/0 . 0/ 

open(l,file=' /dimen dat') 
open(2,file^'/worm dat') 
open (3 , file=' /worml dat') 
open(4, file=Vg ear dat') 
open (5 , file=' /contc dat') 
open (7 , file=' /normal dat') 
open(8 , file=' /pitch dat') 
open(9 , f ile=' /load dat') 

p i=4 . d0*datan ( 1 . dO ) 
bt=bta*pi/180 . dO 
tb=dtan(bt) 
dt=dt2*pi/180 . dO 
df-df l*p i/180. dO 
rfnr=rfn*p i/180 . dO 


calculation of the dimensions 


call dimen( aO , a2 , ratio2 , nl , np , dl , rfnr , bt , 

6c bl , b2 , cc , ddd , dlt , f 20 , f 21 , 

6c rb , rl , r2 , ral , ra2 , rf 1 , rf 2 , rmax2 , 

6c reps , tt2 , cmax, cmin, start , end) 

print*,' DIMEN done' 
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c 

c 

c 

c 

c 


calculation of the worm surfaces along the contact line | 


call worm (a0,a2, rat io,rb,bt,tb, 

& f20, f21 ,bl ,ra2 , rf2 , cc , df , dz) 


c 

c 

c 

c 

c 


c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 


c 

c 

c 

c 


call worml(a0,a2,ratio,rb,bt,tb, 

& f20 , f21 ,bl , ra2 , rf2 , cc , df , dz , reps) 

Print*, 'WORM and WORML done' 


| calcukation of the gear surfaces along contact lines 


t20=f20 

t21=f21 

call gear ( aO , a2 , ratio , rat io2 , rb , tb , bt , rmax2 , cmax , cmin , 
& ra2 , rf 2 , r2 , ral , rf 1 , rl , cc , ddd ,bl , b2 , 

& f 20 , f 21 , df , t20 , t21 , dt , dlt , dz ) 

print*, 'GEAR done' 


calculation of the contact line at the given position tn2 


call contact (aO , a2 , ratio , ratio2 , np , rb , tb , bt , dh , ddd , 

& rl , ral , rf 1 , r2 , ra2 , rf 2 , rmax2 , cmax , cmin , bl , b2 , 

6c f 20 , f21 , df , t20 , t21 , tt2 , tn2 , start , end , dlt , dz , 

6c contc , ncc , nn , exaa , exbb , excc , exdd) 

print* , ' CONTACT done ' 


calculate the characters of pitch points on the worm surface 
at the given position 'tn2' 


call pitch(a0 , a2 , ratio , ratio2 , np , rb , tb ,bt , cc , ddd, 

6c rl , ral , rf 1 , r2 , ra2 , rf 2 , rmax2 , cmax , cmin , bl , b2 , 

6c f 20 , f 21 , df , t20 , t21 , tt2 , tn2 , start , end , dlt , dz , 

6c nn , exaa , exbb , excc , exdd , ep , epp , 

6c tp2,fp2,pxj , pzj ,xpn,ypn, zpn, vx,vy , vz ,ncp) 

print* ,' PITCH done' 


calculate the load distribution at the given position tn2 


if (ncc .ne .ncp) then 

print*, ' ncc. ne. ncp check parameters' 
print* , ncc , ncp 
endif 

call load (a2 , rat io2 , r2 ,np , power , rpm,uf, ncc , tp2 ,pxj , 

6c pzj , xpn , ypn , zpn , contc , wtor , wtor f , f orn , f ornf ) 

print*, 'LOAD done' 

stop 

end 


***************** END of MAIN PROGRAM ****************************** 


c 

c 

c 
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Subroutine DIMEN called by MAIN 

This subroutine is used to develop the dimensions of a wormgear 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 


To run this subroutine the input parameters are: 

aO center distance of the first enveloping (mm) 

a2 center distance of wormgear (mm) 

bta apex angle of the generating plane (deg) 

dl pitch diameter of worm (mm) 

nl number of worm thread 
nbao number of meshing teeth 
ratio speed ratio of wormgear 

rfn normal presure angle (deg) 


The output of the subroutine are: 

bl face width of worm (mm) 

b2 face width of gear (mm) 

cmax maximum throad radius of the gear hob (mm) 

cmin minimum root radius of the gear hob (mm) 

dh clearance (mm) 

dlt half of gear face angle (rad) 

dlta the change in center distance (mm) 

end the angle of the meshing end (rad) 

f20 the angle of the first enveloping begining (rad) 

f21 the angle of the first enveloping ended (rad) 

rl pitch radius of worm (mm) 

r2 pitch radius of gear (mm) 

ral throat radius of worm (mm) 

ra2 throat radius of gear (nun) 

rb radius of base circle (mm) 

rfl root radius of worm (mm) 

rf2 root radius of gear (mm) 

reps the angle of grinding wheel turning to another side (rad) 
rgmax out -radius of the gear (mm) 

start the angle of the meshing starting (rad) 

tao angle between two neighboring gear teeth (rad) 


The results which are not called by other subroutines are written 
in the data file under the name of ' dimen.dat'. 


subroutine dimen ( aO , a2 , ratio , nl , nbao , dl , rfn , b ta , 
& bl,b2,dh,dlta,dlt,f20,f21, 

& rb , rl , r2 , ral , ra2 , rfl , rf 2 , rgmax , 

& reps , tao , cmax , cmin , start , end) 

implicit double precision(a-h,o-z) 

1 format(2x, 6fl2 .4) 

pitch diameter and radii 


n2=idint (nl*ratio) 
d2=2.0*a2-dl 
rl=dl/2 . 0 
r2=d2/2 . 0 


circle pitch, thickness, addendum, dedendum and clearance 


c 

c 

c 
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p i=4 . dO*datan ( 1 . dO ) 
px=pi*d2/n2 
s2=0.55*px 
sl=px- s2 
xm=px/pi 
c 

gama=datan(px*nl/(pi*dl) ) 
r f x=d t an (rfn)/dcos( gama ) 
rfx=datan(rfx) 
dr f =px/ ( 2*d2 ) 
c 

pn=px*dco s ( gama ) 
ba=.225*pn 
bk=0 . 45*pn 
bt=pn/2 . 0 
bf=bt-ba 
dh=bt-bk 
c 

c diameters and radius 
c 

db=d2*ds in ( r f x+dr f ) 
dal=dl+2 . 0*ba 
da2=d2+2 . 0*ba 
dfl=dl+2*ba-2*bt 
df2=d2+2*ba-2*bt 
rb=db/2 . 0 
ral=dal/2 . 0 
ra2=da2/2 . 0 
rfl=dfl/2 . 0 
rf 2=df 2/2 . 0 
c 

rf i=rfx+drf 
tao«=2*pi/n2 

psia=. 5*tao*(ribao~ .45) 
start=rfi-psia 
end=s tar t+2*p s ia 
bl=idint(d2*dsin(psia) ) 
b2=idint( . 9*dfl) 

r gmax=a2 - dsqr t ( rl*r 1 -b2*b2/4 . 0 ) 

cmax=a2-dsqrt(rf2*rf2-bl*bl/4.0) 

cmin=a2 - dsqrt ( ra2*ra2 -bl*bl/4 . 0) 

dlta=cmax+rf 2 - a2 

reps=(2 . 0*rfi+2*sl/d2) 

f20=start-tao 

f21=end+tao 

dl t=das in ( df 1/dl ) 

c : 
c table of the wormgear dimensions 

c 

write(l,*) "Dimensions of selected wormgear" 

write (1,*) "a0,a2,rfn, ratio, db" 

write (1,1) aO, a2,rfn, ratio ,db 

write(l ,*) "p2,s2,sl,bk,bt,cc" 

write(l,l) px,s2,sl,bk,bt,dh 

write ( 1 , *) " dal , df 1 , da2 , df 2 , dl , d2 " 

write (1,1) dal , df 1 , da2 , df 2 , dl , d2 

write (1 , *) " ral , rf 1 , ra2 , rf 2 , rl , r2 ' 

write (1,1) ral,rfl,ra2,rf2,rl,r2 

write ( 1 , *) ' rf i , rfx , drf , reps , lamda , ps ia ' 

write (1,1) rf i , rfx , drf , reps , gama , ps ia 

write (1 , *) 'bta , f 20 , f 21 , start , end' 
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write (1,1) bta , f 20 , f 21 , start , end 

write (1,*) ' dlt,bl,b2' 

write(l,l) dlt,bl,b2 

return 

end 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 


End of Subroutine DIMEN 


Subroutine WORM called by MAIN 

This subroutine is used to calculate the worm surface in SI, where 
worm is ground by a disc . To drive this subroutine , the subroutine 
DIMEN must be called first. 


To run this subroutine the input parameters are: 

aO center distance of the first enveloping (mm) 

a2 center distance of wormgear (mm) 

bl face width of worm (mm) 

bt apex angle of the generating plane (rad) 

cc clearance (mm) 

df step size of fl 

dz step size of zl 

f20 the angle of the first enveloping begining (rad) 

f21 the angle of the second enveloping ended (rad) 

ra2 throat radius of gear (mm) 

ratio speed ratio of wormgear 

rb radius of base circle (mm) 

rf2 root radius of gear (mm) 

tb tan(bt) 


The variables in WORM are: 

fl the running parameter in the first enveloping (rad) 
f2 fl/ratio 

zl the coordinate in worm axis direction 

The results which are not called by other subroutines are written 
in the data file under the name of 'worm.dat'. 


subroutine worm(a0 , a2 , ratio , rb ,bt , tb , 

& f20,f21,bl,ra2,rf2,cc,df ,dz) 

implicit double precision(a-h, o-z) 

2 format (2x, 6fl0 .4) 

Title of W0RM.dat file 

write (2,*) ' The contact lines on the worm surfaces in SI' 

write(2 ,*) ' fl I f2 | zl | xl | yl | 

& curvel' 

write(2,*) ' - 

& ' 


t 


c 

c precalculation 
c 


fi=f20*ratio 
fe=f21*ratio 
m=idnint( (fe-fi)/df)+l 
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zi=*-bl/2 . 0 
l=idint (bl/dz ) 

calculation of the worm surface coordinate in SI 

do 20 k— 1,1 
z=“Z i+(k-l)*dz 
do 10 j-l,m 

fl~fi+(j-l)*df 

f2=fl/ratio 

call eqw( f 1 , f 2 , aO , ratio , rb , tb , 

& p , q , r , s , pi , ql , rl , si , p2 , q2 , r2 , s2 ) 

<iw=p* (q2*rl - ql*r 2 ) +q* (pl*r 2 -p2*rl) +r* (p2*ql -pl*q2 ) 
if(dw.eq.O) goto 10 

fz=- (p*(q2*sl-ql*s2)+q*(pl*s2-p2*sl)+s*(p2*ql-pl*q2) )/dw 
f x= - (q*(r2*sl-s2*rl)+r*(ql*s2-q2*sl)+s*(q2*rl-r2*ql) )/dw 
fy=(p*(r2*sl-s2*rl)+r*(pl*s2-p2*sl)+s*(p2*rl-r2*pl) )/dw 
redge=dsqr t ( f x*fx+f y*fy ) 

dw=pl*q-p*ql 

x= ( ( ql*r - q*r 1 ) *z - q*s l+ql*s ) /dw 
y=_ ( (pl*r-rl*p)*z-p*sl+pl*s)/dw 

judge if the point is in the tooth domain 

if(z.lt.fz) goto 10 

rw=ds qr t ( x*x+y*y ) 

rmax=a2 - dsqr t ( ( r f 2+cc ) **2 - z*z ) 

rmin=a2 - dsqr t ( ( r a2+cc ) **2 - z*z ) 

if (rw. gt.rmax.or.rw.lt. rmin) goto 10, 

justify the tooth undercut 

if (redge.gt.rw) then 

print*, 'the worm surface intersects with the edge curve' 
print*, x,y,z 
endif 

calculation of the worm surface curvature 

f f =p 2*x+q 2*y+r 2*z+s 2 
rr=r*r 

pel=l . 0+p*p/rr 

pfl=p*q/rr 

pgl=l . 0+q*q/rr 

dw=f f *r r*ds qr t ( p*p+q*q+r r ) 

aa=p*rl-r*pl 

bb=q*rl-r*ql 

pll=aa*aa/dw 

pml=aa*bb /dw 

pnl=bb*bb/dw 

curvel=(pel*pnl-2 . 0*pfl*pml+pgl*pll)/ (pel*pgl-pfl*pfl) 

write(2 , 2) fl , f2 , z ,x,y ,curvel 
10 continue 
20 continue 
return 
end 


End of Subroutine WORM 



c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


Subroutine EQW called by WORM 

This subroutine is used to calculate the coefficients of the worm 
surface equations . 


To run this subroutine the input parameters are: 

aO center distance of the first enveloping (mm) 

bt apex angle of the generating plane (rad) 

fl running parameter in the first enveloping (rad) 

f2 fl/ratio 

rb radius of base circle (mm) 

tb tan(bt) 


The output in EQW are: 

p,q,r,s coefficients of the first equation in Eqs(19) 

pl,ql,rl,sl coefficients of the second equation in Eqs(19) 
p2,q2,r2,s2 coefficients of the third equation in Eqs(19) 


subroutine eqw(fl , f 2 , aO , ratio , rb , tb , 

& p , q , r , s , pi , ql , rl , si , p2 , q2 , r2 , s2 ) 

implicit double precis ion(a~h, o-z) 

cf2=dcos(f2) 

sf2=dsin(f2) 

cfl=dcos(fl) 

sfl=dsin(fl) 

Here is the equation (18) 

p=cfl*sf2+tb*sfl 

q=tb*cfl-sfl*sf2 

r=-cf2 

s=rb-a0*sf2 

pl=tb*cfl*ratio~sfl*sf2*ratio+cfl*cf2 
ql=- tb*sf l*ratio - sf l*cf 2 - cf l*sf 2*ratio 
rl=sf2 
sl=-a0*cf2 

p2=-tb*sfl*ratio*ratio~cfl*sf2*ratio*ratio 
& -cfl*sf2-2*sfl*cf2*ratio 

q2=~tb*cfl*ratio*ratio+sfl*sf2*ratio*ratio 
& +sfl*sf2-2*cfl*cf2*ratio 

r2=cf2 
s2=a0*sf2 
return 
end 


=END of Subroutine EQW 


Subroutine WORML called by MAIN 

This subroutine is used to calculate the worm surface in SI, where 
worm is ground by a disc. To drive this subroutine, the subroutine 
DIMEN must be called first. 

To run this subroutine the input parameters are: 


108 



c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 


c 

c 

c 


c 

c 

c 


c 


aO center distance of the first enveloping (mm) 

a2 center distance of wormgear (mm) 

bl face width of worm (mm) 

bt apex angle of the generating plane (rad) 

cc clearance (mm) 

df step size of fl 

dz step size of zl 

f20 the angle of the first enveloping begining (rad) 

f21 the angle of the second enveloping ended (rad) 

ra2 throat radius of gear (mm) 

ratio speed ratio of wormgear 

rb radius of base circle (mm) 

reps angle of the other side of the worm surface (rad) 

rf2 root radius of gear (mm) 

tb tan(bt) 


The variables in WORML are: 

fl the running parameter in the first enveloping (rad) 
f2 fl/ratio 

zl the coordinate in worm axis direction 

The other results which are not called by other subroutines are 
written in the data file under the name of 'worml.dat' 


subroutine worml (aO , a2 , ratio , rb ,bt , tb , 

& f 20 , f 21 ,bl , ra2 , rf2 , cc , df , dz , reps) 

implicit double precision(a-h, o-z) 

3 format(2x, 6fl0.4) 


Title of worml.dat file 


write (3,*) ' 
write(3,*) ' 
write (3,*) ' 
& curve 1' 
write(3,*) ' 
& ' 


The contact lines on the other side worm surfaces' 
fl | f2 j zl { xl ] ylj 


precalculation 

fi=f20*ratio 
fe=f21*ratio 
m=idint( (fe-fi)/df)+l 
zi=-bl/2 . 0 
l=idint (bl/dz ) 

calculation of the worm surface coordinate in SI 


do 20 i=l , 1 
z=zi+(i-l)*dz 
do 10 j=l ,m 

fl=fi+(j-l)*df 

f2=fl/ratio 

f2=f2-reps 

call eqwl ( f 1 , f 2 , aO , ratio , rb , tb , 
i p , q , r , s , pi , ql , rl , si , p2 , q2 , r2 , s2) 

dw=p* (q2*rl - ql*r2)+q* (pl*r2 -p2*rl )+r* (p2*ql -pl*q2 ) 
if (dw.eq.O) goto 10 
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c 


fz— ( p *(q2*sl-ql*s2)+q*(pl*s2-p2*sl)+s*(p2*ql-pl*q2))/dw 

fx~ _ ( q * ( r 2*sl - s2*rl ) +r* ( ql*s2 - q2*sl ) +s* ( q2*rl - r2*ql ) ) /dw 

fy=(p*( r 2 *sl-s 2 *rl)+r*(pl*s 2 -p 2 *sl)+s*(p 2 *rl-r 2 *pl))/dw 

redge=dsqrt(fx*fx+fy*fy) 


c 

c judge if the point is in the tooth domain 
c 

if(z.gt.fz) goto 10 
dw«pl*q-p*ql 

x=( (ql*r-q*rl)*z-q*sl+ql*s)/dw 
y=- ((pl*r-rl*p)*z-p*sl4pl*s)/dw 
rw^dsqrt (x*x+y*y ) 
rmax=a2-dsqrt( (rf2+cc)**2-z*z) 
rmin=a2 - ds qr t ( ( ra2+cc ) **2 - z*z ) 
if (rw. gt.rmax.or.rw.lt. rmin) goto 10 
c 

c judge the tooth undercut 
c 

if (redge.gt.rw) then 

print*, 'this point of worm surface is on the edge curve' 
print*, x,y,z 
endif 
c 

c calculation of the worm surface curverature 
c 

ff— p2*x+q2*y+r2*z+s2 
rr=r*r 

pel=l . 0+p*p/rr 

pfl=p*q/rr 

pgl=l . 0+q*q/rr 

dw=ff*rr*dsqrt(p*p+q*q+rr) 

aa=p*rl-r*pl 

bb=q*rl-r*ql 

pll=“aa*aa/dw 

pml=aa*bb / dw 

pnl=bb*bb /dw 

curvel=(pel*pnl-2 .0*pfl*pml+pgl*pll)/(pel*pgl-pfl*pfl) 
c 

write (3,3) f 1 , f 2+reps , z , x , y , curvel 
10 continue 

20 continue 
return 
end 
c 

c ■ ■ = — End of Subroutine W0RML — 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


Subroutine EQWL called by W0RML 

This subroutine is used to calculate the coefficients of the worm 
surface . 


To run this subroutine the input parameters are: 

aO center distance of the first enveloping (mm) 

bt apex angle of the generating plane (rad) 

fl running parameter of the first enveloping (rad) 

f2 fl/ratio 

rb radius of base circle (mm) 

tb tan(bt) 
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c 

c The output in EQW are: 
c 

c p,q,r,s coefficients of the worm equation 

c pl,ql,rl,sl coefficients of the worm equation 

c p2,q2,r2,s2 coefficients of the worm equation 

c 

c 

c 

subroutine eqwl ( f 1 , f 2 , aO , ratio , rb , tb , 

& p , q , r , s , pi , ql , rl , si , p2 , q2 , r2 , s2) 

implicit double precision(a-h , o-z) 

cf2=dcos(f2) 

sf2=dsin(f2) 

cfl=dcos(fl) 

sfl=dsin(fl) 

c 

p=cfl*sf2+tb*sfl 
q=tb*cfl-sfl*sf2 
r=-cf 2 

s=-rb-a0*sf2 

pl=tb*cf l*ratio - sf l*sf 2*ratio+cf l*cf 2 

ql=-tb*sfl*ratio-sfl*cf2-cfl*sf2*ratio 

rl=sf2 

sl=-a0*cf2 

p2=- tb*sf l*ratio*ratio - cf l*sf 2*ratio*ratio 
& -cfl*sf2-2*sfl*cf2*ratio 

q2=-tb*cfl*ratio*ratio+sfl*sf2*ratio*ratio 
& +sfl*sf2-2*cfl*cf2*ratio 

r2=cf2 
s2=a0*sf2 
return 
end 


c 




c = 

c 

c 


• " Cilia. OX D UDi OULlIic CiV^Wi-i 


c = 
c 

Subroutine 

GEAR called by MAIN 


c 

This subroutine is used to calculate the gear surface 

in S2 , where 

c 

the gear hob is ground by a disc. To drive this subroutine, the 

c 

subroutine 

DIMEN must be called first. 


c 

c 

To run this subroutine the input parameters are: 


c 

aO 

center distance of the first enveloping 

(mm) 

c 

a2 

center distance of wormgear 

(mm) 

c 

bl 

face width of worm 

(mm) 

c 

b2 

face width of gear 

(mm) 

c 

bt 

apex angle of the generating plane 

(rad) 

c 

cc 

clearance 

(mm) 

c 

cmax 

maximum throad radius of the gear hob 

(mm) 

c 

cmin 

minimum root radius of the gear hob 

(mm) 

c 

ddd 

the change in center distance 

(mm) 

c 

dlt 

half of gear face angle 

(rad) 

c 

df 

step size of fl 


c 

dz 

step size of z2 


c 

f20 

the angle of the first enveloping begining 

(rad) 

c 

f21 

the angle of the second enveloping end 

(rad) 

c 

rl 

pitch radius of worm 

(mm) 

c 

r2 

pitch radius of gear 

(mm) 

c 

ral 

throat radius of worm 

(mm) 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 


c 


c 

c 

c 


c 

c 

c 


c 

c 

c 


ra2 throat radius of gear (mm) 

ratio speed ratio of wormgear 

ratio2 speed ratio of wormgear 

rb radius of base circle (mm) 

rfl root radius of worm (mm) 

rf2 root radius of gear (mm) 

rmax2 out -radius of the gear (mm) 


t20 running angle of the second enveloping start (rad) 

t21 running angle of the second enveloping ended (rad) 

tb tan(bt) 

The variables in GEAR are: 

tl the running parameter in the second enveloping (rad) 
t2 tl/ratio 

z2 the coordinate in gear axis direction 

The results which are not called by other subroutines are written 
in the data file under the name of 'gear.dat'. 


subroutine gear (aO, a2 , ratio ,ratio2 , rb , tb ,bt , rmax2 , cmax, cmin, 

& ra2 , r f 2 , r 2 , ral , rfl , r 1 , cc , ddd , b 1 , b2 , 

& f 20 , f 21 , df , t20 , t21 , dt , dlt , dz ) 

implicit double precision(a-h, o-z) 

4 format(2x, 6fl2 .4) 

Title of Gear.dat file 

write(4,*)' This is the point coordinates on gear surface: ( when 
Sc tl— fl, the point is on the first contact line)' 
write (4,*) ' 

write (4,*) ' tl | fl | x2 | y2 | 

& z2 I r' 

write (4,*) ' 

& - ' 

l=idint (b2/dz ) +1 
m=idint( (t21- t20)/dt)+l 
dfl=df/10.0 

n=idint(4 . dO*dlt/dfl)+l 

precalculation 

do 10 i=l,m 

t2=t20+(i-l)*dt 

tl=t2*ratio2 

calculation of the curved part of the gear surface in S2 
do 20 j=l,n 

fl=(tl-2.0*dlt)+(j-l)*dfl 

f2=fl/ratio 

if (f2 . It . f20 . or . f2 . gt . f21) goto 20 

call second(t2 , f2 , aO , a2 , ratio , ratio2 , rb , tb ,bt , 

Sc xl ,yl , zl ,x2 ,y2 , z2 ,xn,yn, zn, 

& curve! , curve 2 , curve , dw) 

judge if the point is in the tooth domain 
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if(dw.eq.O) goto 20 

if (curve. It. 0) goto 20 

if(z2.gt.b2/2. or. z2.lt. -b2/2) goto 20 

rg=dsqrt(x2*x2+y2*y2) 

rmin== ( a2+ddd) - dsqr t ( cmax**2 - z2*z2 ) 

rmax= ( a2+ddd) - ds qr t ( cmin**2 - z 2*z 2 ) 

if ( rmax . gt . rmax2 ) then 

rmax=rmax2 

endif 

if (rg. gt. rmax. or .rg.lt. rmin) goto 20 

write (4 , 4) tl , f 1 , x2 , y2 , z2 , rg 
20 continue 


c 

c 

c 

c 


c 

c 

c 


c 


c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


if (a0.ne.a2.or.ratio.ne.ratio2) goto 10 

calculation of the planar part of the gear surface in S2 
( for the case of fl=tl in the conventional form ) 

do 30 k=l , 1 

z2=-b2/2+(k-l)*dz 

fl=tl 

f2=t2 

if(f2.1t.f20.or.f2.gt.f21) goto 30 

call f irst ( t2 , f2 , aO , a2 , ratio , ratio2 , rb , tb ,bt , z2 , 

& xl,yl,zl,y2,x2,xn,yn,zn, 

& curve 1 , curve 2 , curve , dw ) 

judge if the point is in the tooth domain 

if(dw.eq.O) goto 30 

if (curve. It. 0) goto 30 

rg=dsqrt (x2*x2+y2*y2 ) 

rmin= ( a2+ddd) - dsqr t ( cmax**2 - z2*z 2 ) 

rmax= ( a2+ddd) - dsqr t ( cmin**2 - z 2*z 2 ) 

if ( rmax . g t . rmax2 ) then 

rmax=rmax2 

endif 

if (rg. gt.rmax.or.rg.lt. rmin) goto 30 

write (4 , 4) tl , f 1 , x2 , y2 , z2 , rg 
30 continue 

10 continue 
return 
end 


End of Subroutine GEAR 


Subroutine SECOND called by Subroutine GEAR 

This subroutine is used to calculate the coefficiences of the gear 
surface equations Eq(28) in techanical report. 


To run this subroutine the input parameters are: 

aO center distance of the first enveloping (mm) 
a2 center distance of wormgear (mm) 
btr apex angle of the generating plane (rad) 
f2 running parameter in the first enveloping (rad) 
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c 

ratio 

speed ratio of wormgear 

c 

ratio2 

speed ratio of wormgear 

c 

rb 

radius of base circle 

c 

t2 

running parameter in the second 

c 

c 

tb 

tan(bt) 

c 

c 

The output 

in 

SECOND are: 

c 

c 

curvel 


worm surface curvature at 

c 

curve 2 


gear surface curvature at 

c 

curve 


relative curvature at the 

c 

dw 


index ( dw=0 the point is 

c 

xl ,yl , zl 

contact point coordinates 

c 

x2,y2,z2 

contact point coordinates 

c 

xn , yn , 

zn 

contact point coordinates 

c • 
c 





(mm) 

(rad) 


, *■«- j «-v , w**- y , - 7 7 t ' 

xl , yl , zl , x2 , y 2 , z2 , xn , yn , zn , 

curve! , curve 2 , curve , dw) 
implicit double precision(a-h, o-z) 

calculate the coordinates of the point for given fl, tl 

tl=t2*ratio2 

fl=f2*ratio 

call eqg ( tl , t2 , f 1 , f 2 , aO , a2 , ratio , ratio2 , tb , rb , b tr , 

& p,q,r,s,pl,ql,rl,sl,p2,q2,r2,s2) 

dw«p*(q2*rl-ql*r2)+q*(pl*r2-p2*rl)+r*(p2*ql-pl*q2) 

if (dw.eq.O) goto 10 

x2=- (q*(r 2 *sl-s 2 *rl)+r*(ql*s 2 -q 2 *sl)+s*(q 2 *rl-r 2 *ql))/dw 
y 2 =(p*(r 2 *sl-s 2 *rl)+r*(pl*s 2 -p 2 *sl)+s*(p 2 *rl-r 2 *pl))/dw 
z2—~ (p*(q2*sl-ql*s2)+q , *(pl*s2-p2*sl)+s'*(p2*ql-pl*q2) )/dw 

ctl=dcos (tl) 

stl=dsin(tl) 

ct2=dcos(t2) 

st2=ds in(t2) 

cf2=dcos(f2) 

sf2=dsin(f2) 

cfl=dcos(fl) 

sfl=dsin(fl) 

sft=dsin(fl- tl) 

cft=dcos(fl-tl) 

Coordinate transformation 

xl=-x2*ctl*ct2+y2*ctl*s t2 - z2*s tl+a2*ctl 

yl=x2*stl*ct2-y2*stl*st2-z2*ctl-a2*stl 

zl=-x2*st2 -y2*ct2 

xn=x2*ct2 -y2*st2 

yn=x2*st2+y2*ct2 

zn=z2 

curvature of the gear surface 


rr=r*r 

pe2=l+p*p/rr 

pf2=p*q/rr 

pg2=l+q*q/rr 
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c 


b=*p*rl - r*pl 
d=q*rl-r*ql 
c=p*r2-r*p2 
e=q*r2-r*q2 

rto=ratio*ratio 

rto2=ratio2*ratio2 

f fx=sf 2*ct2*cf t+tb*ct2*sf t+2 . 0*cf2*ct2*sf t/ratio 
& +(sf2*ct2*cft-cf2*st2)/rto 

f fy-sf 2*st2*cf t+tb*st2*sf t+2 . 0*cf 2*st2*sf t/ratio 
& +(sf2*st2*cft+cf2*ct2)/rto 

ffz=sf 2 *sft-tb*cft- 2 . 0 *cf 2 *cft/ratio+sf 2 *sft/rto 
ffd=a2*(sf2*cft+tb*sft+2 . 0*cf2*sft/ratio+sf2*cft/rto) 

6c -a0*sf2/rto 

f f =f fx*x2 - f fy*y2 - f f z*z2 - f f d 

ftx=sf 2 *st 2 *sft-tb*st 2 *cft+ratio 2 *(sf 2 *ct 2 *cft+tb*ct 2 *sft) 

6c +(sf 2 *ct 2 -cf 2 *st 2 *cft+ratio 2 *cf 2 *ct 2 *sft)/ratio 

f ty =-sf 2 *ct 2 *sft+tb*ct 2 *cft+ratio 2 *(sf 2 *st 2 *cft+tb*st 2 *sft) 
6c +(sf 2 *st 2 +cf 2 *ct 2 *cft+ratio 2 *cf 2 *st 2 *sft) /ratio 
ftz-ratio2*(sf2*sft-tb*cft-cf2*cf t/ratio) 
ftd=ratio2*a2*(sf2*cft+tb*sft+cf2*sf t/ratio) 
f t=- f tx*x2+f ty*y2+f tz*z2+f td 
c 

ttx=-p+2.0*ratio2*(sf2*st2*sft-tb*st2*cft) 

6c +rto2*(sf2*ct2*cft+tb*ct2*sft) 

tty=-q+2 . 0*ratio2* ( sf 2*ct2*sf t - tb*ct2*cf t ) 

6c -rto2*(sf2*st2*cft+tb*st2*sft) 

ttz=rto2* ( tb*cf t - s f 2*sf t ) 
ttd=rto2*a2*(sf2*cft+tb*sft) 
tt=ttx*x2+tty*y2+ttz*z2 - ttd 
c 

g=ff*tt-ft*ft 

cm=dcos(btr)/(g*rr) 

pl2=cm* (b*b*tt - 2*b*c*f t+c*c*f f ) 

pm2=cm* (b*d*tt - (b*e+c*d) *f t+c*e*f f ) 

pn2=cm* ( d*d* tt - 2*d*e*f t+e*e*f f ) 

curve2= (pe2*pn2 - 2*pf 2*pm2+pg2*pl2 ) / (pe2*pg2 -pf 2*pf 2) 
c 

c curvature of the worm surface 
c 

pel=l . 0+p*p/rr 

pfl=p*q/rr 

pgl=l . 0+q*q/rr 

dw=f f *rr*dsqrt (p*p+q*q+rr ) 

aa=p*rl-r*pl 

bb=q*rl-r*ql 

pll=aa*aa/dw 

pml=aa*bb /dw 

pnl=bb*bb/dw 

curvel=(pel*pnl-2.0*pfl*pml+pgl*pll)/(pel*pgl-pfl*pfl) 

c 

c relative curvature 
c 

curve=curve 2 - curve 1 
c 

10 return 
end 


End of Subroutine SECOND 


c 

c 

c 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


c 

c 

c 


c 


Subroutine EQG called by Subroutine FIRST and SECOND 

This subroutine is used to calculate the coefficients of the gear 
surface equations. 


To run this subroutine the input parameters are: 

aO center distance of the first enveloping (mm) 

a2 center distance of the first enveloping (mm) 

bt apex angle of the generating plane (rad) 

fl running parameter of the first enveloping (rad) 

f2 fl/ratio 

tl running parameter of the second enveloping (rad) 

t2 tl/ratio 

ratio speed ratio of wormgear 
ratio2 speed ratio of wormgear 

rb radius of base circle (mm) 

tb tan(bt) 


The output in EQW are : 

p,q,r,s coefficients in Eqs(29) 
pl,ql,rl,sl coefficients in Eqs(30) 
p2,q2,r2,s2 coefficients in Eqs(31) 


subroutine eqg ( tl , t2 , f 1 , f 2 , aO , a2 , ratio , rat io2 , tb , rb , btr , 

& p , q , r , s , pi , ql , rl , si , p2 , q2 , r2 , s2 ) 

implicit double precision(a-h.o-z) 

ctl=dcos (tl) 

stl=dsin(tl) 

ct2=dcos(t2) 

st2=dsin(t2) 

cf2=dcos(f2) 

sf2=dsin(f2) 

cfl=dcos(fl) 

sfl=dsin(fl) 

sft— dsin(fl- tl) 

cft=dcos(fl-tl) 

Here is the Eqs(29)-(30) 

p=-sf2*ct2*cft-tb*ct2*sft+cf2*st2 
q=sf2*st2*cft+tb*st2*sft+cf2*ct2 
r=sf2*sft-tb*cft 
s =a2*(sf2*cft+tb*sft) -a0*sf2+rb 

pl=sf 2 *ct 2 *sft-tb*ct 2 *cft-(sf 2 *st 2 +cf 2 *ct 2 *cft)/ratio 
ql=-sf2*st2*sft+tb*st2*cft - ( sf 2*ct2 - cf 2*s t2*cft) /ratio 
rl=sf2*cft+tb*sft+cf2*sft/ratio 

sl=-a2*(sf2*sft-tb'*cft-cf2*cft/ratio) -a0*cf2/ratio 
p2=q - ratio2* (sf 2*c t2*sf t - tb*ct2*cf t ) 
q2=-p+ratio2*(sf2*st2-*sft’tb*st2*cft) 
r2= - ratio2* ( s f 2*cf t+tb*sf t ) 
s2=a2*ratio2*(sf2*sft-tb*cft) 

return 

end 


End of Subroutine EQG 


c 

c 

c 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 


c 


c 

c 

c 


Subroutine FIRST called by GEAR 

This subroutine is used to calculate the coefficients of the gear 
surface equations and curvatures. 


To run this subroutine the input parameters are: 

aO center distance of the first enveloping (mm) 

a2 center distance of wormgear (nun) 

btr apex angle of the generating plane ^ (tad) 

f2 running parameter in the first enveloping (rad) 

ratio speed ratio of wormgear 
ratio2 speed ratio of wormgear 

rb radius of base circle (mm) 

t2 running parameter in the second enveloping (rad) 

tb tan(bt) 

z2 coordinate of the point in S2 


The output in FIRST are : 


curve 1 
curve2 
curve 
dw 

xl ,yl ,zl 
x2 ,y2 
xn,yn,zn 


worm surface curvature at 
gear surface curvature at 
relative curvature at the 
index ( dw=0 the point is 
contact point coordinates 
contact point coordinates 
contact point coordinates 


the contact point 
the contact point 
contact point 
a singularity point ) 
in SI 
in S2 
in Sn 


subroutine f irs t ( t2 , f 2 , aO , a2 , ratio , ratio2 , rb , tb , btr , z2 , 
& xl,yl,zl,y2,x2,xn,yn,zn, 

& curvel , curve2 , curve , dw) 

implicit double precision(a-h, o-z) 


calculation of the coordinates of the point 
at the given position fl, tl 

tl=t2*ratio2 

fl=f2*ratio 

call eqg ( tl , t2 , f 1 , f 2 , aO , a2 , ratio , rat io2 , tb , rb , b tr , 

& p , q , r , s , pi , ql , rl , si , p2 , q2 , r 2 , s2 ) 

dw=pl*q - p*ql 
if (dw.eq.O) goto 10 
x2= ( ( ql*r - q*rl ) *z 2 - q*s l+ql*s ) /dw 
y2=- ( (pl*r-p*rl)*z2-p*sl+pl*s)/dw 


ctl=dcos(tl) 

stl=dsin(tl) 

ct2=dcos(t2) 

st2=*dsin(t2) 

cf2=dcos(f2) 

sf2=dsin(f2) 

cfl=dcos(fl) 

sfl=dsin(fl) 

sft=dsin(fl- tl) 

cft=dcos (fl-tl) 


Coordinate transformation 
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xl— x 2 *ctl*ct 2 +y 2 *ctl*st 2 -z 2 *stl+a 2 *ctl 
yl=x2*stl*ct2 -y2*s tl*s t2 - z2*ctl - a2*stl 
zl=-x2*st2-y2*ct2 
xn=x2*ct2 -y2*st2 
yn==x2*st2+y2*ct2 
zn=z2 
c 

c curvature of the gear surface 


rr=r*r 

pe2=l+p*p/rr 
pf2=p*q/rr 
pg2— l+q*q/rr 

b=p*rl-r*pl 

d=q*rl-r*ql 

c=b 

e=d 

rto=ratio*ratio 

rto2=ratio2*ratio2 

f fx=sf 2*ct2*cf t+tb*ct2*sf t+2 . 0*cf 2*ct2*sf t/ratio 
& +(sf2*ct2*cft-cf2*st2)/rto 

f fy»sf 2*st2*cf t+tb*st2*sf t+2 . 0*cf2*st2*sft/ratio 
& + ( sf 2*s t2*cf t+cf 2*ct2 ) /rto 

ffz=sf2*sft-tb*cft-2.0*cf2*cft/ratio+sf2*sft/rto 

ffd=a2*(sf2*cft+tb*sft+2.0*cf2*sft/ratio+sf2*cft/rto) 

& -a0*sf2/rto 

ff=ffx*x2-ffy*y2-ffz*z2-ffd 

ftx=sf2*st2*sft-tb*st2*cft+ratio2*(sf2*ct2*cft+tb*ct2*sft) 

& + ( sf 2*ct2 - cf 2*s t2*cf t+ratio2*cf 2*ct2*sf t) /ratio 

fty=-sf 2 *ct 2 *sft+tb*ct 2 *cft+ratio 2 *(sf 2 *st 2 *cft+tb*st 2 *sft) 
& +(sf 2 *st 2 +cf 2 *ct 2 *cft+ratio 2 *cf 2 *st 2 *sft) /ratio 
ftz=ratio2*(sf2*sft-tb*cft-cf2*cf t/ratio) 
ftd=ratio2*a2* (sf2*cft+tb*sf t+cf 2*sf t/ratio) 
f t=- f tx*x2+f ty*y2+f tz*z2+f td 

ttx=-p+2 . 0*ratio2*(sf2*st2*sft-tb*st2*cft) 

& +rto2*(sf2*ct2*cft+tb*ct2*sft) 

tty==-q+2 . 0*ratio2*(sf2*ct2*sft-tb*ct2*cft) 

& -rto2*(sf2*st2*cft+tb*st2*sft) 

ttz=rto2*(tb*cft-sf2*sft) 
ttd=rto2*a2*(sf2*cft+tb*sft) 
tt=ttx*x2+tty*y2+ttz*z2 - ttd 
c 

g=ff*tt-ft*ft 

cm=dcos(btr)/(g*rr) 

pl2=0 

pm2=0 

pn2=0 

curve2= (pe2*pn2 - 2*pf 2*pm2+pg2*pl2 ) / (pe2*pg2 -pf 2*pf 2 ) 
c 

c curvature of the worm surface 
c 

pel=l . 0+p*p/rr 

pfl=p*q/rr 

pgl=l . 0+q*q/rr 

dw=f f *rr*dsqrt (p*p+q*q+rr ) 

aa=p*rl-r*pl 

bb=q*rl-r*ql 
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c 

c 

c 


pll=aa*aa/dw 
pml=aa*bb /dw 

pnl=bb*bb/dw „ 

curve 1= ( pe l*pnl - 2 . 0*pfl*pml+pgl*pll)/(pel*pgl-pfl*pfl) 

relative curvature 

curve=curve 2 - curve 1 


c 

10 return 
end 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


END of Subroutine FIRST 


Subroutine CONTACT called by MAIN 

This subroutine is used to calculate the real contact lines in Sn 
for the given position tn2. To drive this subroutine, the 
subroutine DIMEN must be called first and the given position must 
be given previously. 


To 


run this subroutine the input parameters are: 
aO center distance of the first enveloping 

a2 center distance of wormgear 

bl face width of worm 

b2 face width of gear 

btr apex angle of the generating plane 
dh clearance 

cmax maximum throad radius of the gear hob 

cmin minimum root radius of the gear hob 

ddd the change in center distance 

dlt half of gear face angle 

df step size of fl 

dz step size of z2 

end the angle of the meshing ended 

f20 the angle of the first enveloping begining 

f21 the angle of the second enveloping ended 

np number of meshing teeth 

rl pitch radius of worm 

r2 pitch radius of gear 

ral throat radius of worm 

ra2 throat radius of gear 

ratio speed ratio of wormgear 

ratio2 speed ratio of wormgear 

rb radius of base circle 

rfl root radius of worm 

rf2 root radius of gear 

rmax2 outradius of the gear 

start the angle of the meshing starting 

t20 running angle of the second enveloping start 

t21 running angle of the second enveloping ended 

tb tan(bt) 

tn2 the given position 

tt2 angle between two neighboring gear teeth 


(mm) 

(mm) 

(mm) 

(mm) 

(rad) 

(mm) 

(mm) 

(mm) 

(mm) 

(rad) 


(rad) 

(rad) 

(rad) 

(mm) 

(mm) 

(mm) 

(mm) 


(mm) 

(mm) 

(mm) 

(mm) 

(rad) 

(rad) 

(rad) 


(rad) 


The output of the subroutine are: 

contc the length of contact lines of each meshing tooth 
nc the number of real meshing pair 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 


c 

c 

c 


c 

c 

c 

c 


The variables in CONTACT are: . . 

tl the running parameter in the second enveloping (rad; 
t2 tl/ratio 

z2 the coordinate in gear axis direction 


The common exchange dimensions are: __ 

x y,z save the coordinates of contact point on meshing pair 
nn maximum dimension number of the exchange dimension 

rg save the radius of contact point to gear axis 


The results which are not called by other subroutines are written 
in the data file under the name of 'contc.dat' . 


subroutine contact (aO , a2 , ratio , ratio2 ,np , rb , tb ,btr , dh, ddd, 

& rl , ral , rf 1 , r2 , ra2 , rf 2 , rmax2 , cmax , cmin,bl ,b2 , 

& f 20 , f 21 , df , t20 , t21 , tt2 , tn2 , start , end , dlt , dz , 

& contc,nc,nn,x,y,z ,rg) 

implicit double precision(a-h,o-z) 
dimension x(nn) ,y(nn) ,z(nn) ,rg(nn) ,contc(np) 

5 format (2x, 6fl0 .4) 

Title of contc.dat file 

write(5,*)' The length of contact lines on each meshing teeth' 

write (5,*) ' r T ” i 

write (5,*) ' tnl | first | second j total 

write(5 ,*) ' ' 

precalculation 

pi=4 . d0*datan(l . dO) 
dfl=df/100 . 0 
dz2=dz/10 . 0 
nc=0 

l=idint (b2/dz 2 ) +1 
m=idint ( f 21/tt2) +1 
n==idint (4 . 0*dlt/dfl)+l 


calculation of the coordinates of the contact point 

do 10 i-l,m 

t2=tn2+ ( i - 1 ) *tt2 

if (t2.lt. start. or .t2.gt. end) goto 10 

tl=ratio2*t2 

ctl=dcos(tl) 

stl=dsin(tl) 

ct2=dcos(t2) 

st2=dsin(t2) 

ind=0 


do 20 i=l ,n 

fl=(tl-2.0*dlt)+(j-l)*dfl 

f2=fl/ratio 

if ( f 2. It. start. or. f 2. gt. end) goto 20 
call second (t2 , f2 , aO , a2 , ratio , ratio2 , rb , tb ,bt , 
& xl,yl,zl,x2,y2,z2,xn,yn,zn, 
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6c curvel , curve2 , curve , dw) 

judge if the point is in the tooth domain 

if(dw.eq.O) goto 20 
if (curve. It. 0) goto 20 

if(z2.1t. -b2/2.0.or.z2.gt.b2/2.0) goto 20 
r g2=dsqr t ( x2*x2+y2*y 2 ) 

rn2==(a2+ddd) -dsqrt( (cmax-dh)**2-z2*z2) 

rm2-= ( a2+ddd) - dsqr t ( cmin**2 - z2*z2 ) 

if ( rm2 . gt . rmax2 ) then 

rm2=rmax2 

endif 

if (rg2. gt.rm2.or.rg2.lt. m2) goto 20 

if (zl . gt .bl/2 . 0 . or . zl . It . -bl/2 . 0) goto 20 

rwl=ds qr t ( xl*xl+y l*y 1 ) 

rnl— a2 - dsqrt ( (ra2+dh) **2 - zl*zl ) 

rml=a2 - dsqrt ( ( rf 2+dh) **2 - zl*zl) 

if (rwl.gt. rml.or.rwl.lt. rnl) goto 20 

ind=ind+l 

if (ind. gt .nn) then 

print *, ' demens ions nn is not big enough for contact' 
print*,' increase nn or df in main program. ' 
print* , ind 
pause 
endif 
x(ind)=xn 
y(ind)=yn 
z (ind)=z2 
rg(ind)=rg2 
20 continue 

calculation of the length of the second contact line 
if (ind.ne.0) then 

rcl2=dsqrt( (x(ind) -x(l) )**2+(y(ind) -y(l) )**2 
& +(z(ind) -z(l) )**2) 

endif 

do 30 j=l,ind 
x(j)=0 

y(j)=° 

z(j)=0 

30 continue 


if (a0.ne.a2 . or. ratio. ne.ratio2) goto 10 

calculation of the coordinates of the contact point 
(this is for the situation where fl=tl) 

ind=0 

do 60 k=l , 1 

z2=-b2/2 . 0+(k-l)*dz2 

fl=tl 

f2=t2 

if(f2.1t.f20.or.f2.gt.f21) goto 60 

call first(t2 , f2 ,a0 , a2 , ratio , ratio2 , rb , tb ,bt , z2 , 
& xl,yl,zl,y2 ,x2 ,xn,yn,zn, 

& curvel , curve2 , curve , dw) 
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c 

c 

c 


judge if the point is in the tooth domain 


c 


c 


c 


c 

c 

c 


c 

c 

c 


c 


c 


if(dw.eq.O) goto 60 
if (curve. It. 0) goto 60 

if(zl.gt.bl/2.0. or. zl.lt. -bl/2.0) goto 60 
rwl=dsqrt (xl*xl+yl*yl ) 
rml«»a2- dsqrt ( (rf2+dh)*(rf2+dh) -zl*zl) 
rnl=a2 - dsqrt ( (ra2+dh)*(ra2+dh) -zl*zl) 
if (rwl. gt.rml.or.rwl.lt. ml) goto 60 

if (z2 . It . -b2/2 . 0 . or . z2 . gt .b2/2 . 0) goto 60 
rg2=dsqrt(x2*x2+y2*y2) 

rn2=(a0+ddd) - dsqrt ( (cmax-dh)**2-z2*z2) 

rm2= ( a0+ddd) - ds qr t ( cmin**2 - z 2*z2 ) 

if (rm2.gt.rmax2) then 

rm2=rmax2 

endif 

if (rg2. gt.rm2.or.rg2.lt. m2) goto 60 
ind=ind+l 

if (ind. gt ,nn) then 

print*, ' demensions nn is not big enough for contact' 
print*,' increase nn or dz in main program. ' 
print* , ind 
pause 
endif 
x(ind)=xn 
y(ind)=yn 
z(ind)=z2 
rg(ind)=rg2 
60 continue 

calculation of the length of the second contact line 
if (ind.ne.0) then 

rcll=dsqrt( (x(ind) -x(l) )**2+(y(ind) -y(l))**2 
& +(z (ind) -z(l) )**2) 

endif 

calculation of the sum of two contact lines on the meshing tooth 

if(rcll.eq.0.and.rcl2.eq.0) goto 10 
nc=nc+l 

contc (nc)=rcll+rcl2 

write (5 , 5) tl,rcll,rcl2, contc (nc) 

rcll=0 

rcl2=0 

do 70 j=l,ind 
x(j)=0 

y(j)=o 

z(j)-0 

70 continue 
10 continue 


return 

end 

-^■= END of Subroutine CONTACT 
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c 

c 

c = 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


Subroutine PITCH called by MAIN 

This subroutine is used to calculate the coordinates of the wora 
pitch points, its normal direction and relative velocity at this 
point in S j . All the calculations are for the given position tn2 
To drive this subroutine DIMEN must be called first and the given 
position must be input previously. 


To run this subroutine the input parameters are: 

aO center distance of the first enveloping 

a2 center distance of wormgear 

bl face width of worm 

b2 face width of gear 

btr apex angle of the generating plane 

dh clearance 

cmax maximum throat radius of the gear hob 

cmin minimum root radius of the gear hob 

ddd the change in center distance 

dlt half of gear face angle 

df step size of fl 

dz step size of z2 

end the angle of the meshing ended 

ep,epp convergence criterion for root searching 

f20 the angle of the first enveloping begining 

f21 the angle of the second enveloping ended 

np number of meshing teeth 

rl pitch radius of worm 

r2 pitch radius of gear 

ral throat radius of worm 

ra2 throat radius of gear 

ratio speed ratio of wormgear 

ratio2 speed ratio of wormgear 

rb radius of base circle 

rfl root radius of worm 

rf2 root radius of gear 

rmax2 outradius of the gear 

start the angle of the meshing starting 

t20 running angle of the second enveloping start 

t21 running angle of the second enveloping ended 

tb tan(bt) 

tn2 the given position 

tt2 angle between two neighboring gear teeth 


(mm) 

(mm) 

(mm) 

(ram) 

(rad) 

(mm) 

(mm) 

(mm) 

(mm) 

(rad) 


(rad) 

(rad) 

(rad) 

(mm) 

(mm) 

(mm) 

(mm) 


(mm) 

(mm) 

(mm) 

(mm) 

(rad) 

(rad) 

(rad) 


(rad) 


The 


output 
nc 
tp2 
fp2 

pxj ,pzj 
xpn , ypn , zpn 

vx,vy ,vz 


of the subroutine are: 

the number of real meshing teeth 
the position of each meshing pair 
the f2 value of pitch points 
the coordinates of pitch points in Sj 
the normal direction components of 
pitch point in Sj 

the relative velocity components of 
pitch point in Sj 


The variables in PITCH are: 

tl the running parameter in the second enveloping (rad) 
t2 tl/ratio 
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zl the coordinate in worm axis direction 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 


The common exchange dimensions are: 

ff save the £1 values of the worm pitch points 

xm,zm save the coordinates of the worm mid-plane profile 

rr save the radius of worm pitch points 

nn maximum dimension number of the exchange dimension 


The results which are not called by other subroutines are written 
in the data file under the name of 'pitch.dat'. 


subroutine pitch(aO , a2 , ratio , rat io2 ,np , rb , tb ,btr , dh, ddd, 

& rl , ral , r f 1 , r 2 , ra2 , rf 2 , rmax2 , cmax , cmin , bl ,b2 , 

& f 20 , f 21 , df , t20 , t21 , tt2 , tn2 , start , end, dlt , dz , 

& nn,xm,zm,ff ,rr,ep,epp, 

& tp2 , fp2 ,pxj ,pzj , xpn , ypn , zpn , vx , vy , vz , nc ) 

implicit double precision(a-h, o-z) 

dimension xm(nn) ,zm(nn) ,ff(nn) ,rr(nn) ,tp2(np) ,fp2(np) , 

& pxj (np) ,pzj (np) ,xpn(np) ,ypn(np) ,zpn(np) , 

& vx(np) ,vy(np) ,vz(np) 

6 format (2x, 6fl2 .4) 


Title of contc . dat file 

write(7,*) ' pitch point and normal direction ' 

write (7 ,*) ' ; 7 7 7 

write (7,*) ' tnl | xj | zj | xpn 

& ypn I zpn' 

write( / ,*) ' - 


write (8 ,*) ' pitch point and relative velocity direction', 

write(8 ,*) ' " “ “ 7 : 7 7 7 7 ~ ’ 

write (8,*) ' tnl | xj [ zj | rvx | 

& rvy I rvz' 

write (8 ,*) ' 


c Precalculation 
c 

nc==0 
yl=0 . 0 

p i=4 . d0*datan ( 1 . dO ) 
m=idint ( f 21/tt2 ) +1 
startl=start*ratio 
endl=end*ratio 
df l=df /100 .0 

ml=idint(pi/ (2 . 0*dfl) )+l 
c 

c Worm profile on the mid-plane (yj=0) 
c 

do 50 k=l,m 
ind=0 

fn2=tn2+ (k- 1 ) *tt2 
fnl=fn2*ratio2 

if ( fnl . It . startl . or . fnl . gt . endl ) goto 50 
c 

do 10 i=l,ml 

f l=f nl -pi/4.0+(i-l)*dfl 
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call wo rat < f 1 , fnl , yl , aO , ratio , rb , tb , b tr , xl , zl , mw) 
if(mw.eq.l) goto 10 
if(xl.lt.O) goto 10 
if (dabs(zl) .gt.rb) goto 10 

rwl=dsqrt(xl*xl+yl*yl) 
rnl=aO-dsqrt( (ra2+dh)*(ra2+dh) -zl*zl) 
rml=a0 - dsqrt ( (rf2+dh)*(rf2+dh) -zl*zl) 
if (rwl. gt.nnl.or.rwl.lt. ml) goto 10 

ind=ind+l 

i f ( ind . gt . nn) then 

print*, 'the dimension is not big enough in PITCH' 

print* ,' increase nn in main' 

print*, ind 

pause 

endif 

xm(ind)=xl 
zm(ind)— zl 

rr ( ind) =dsqr t ( ( aO -xl) * ( aO -xl) +zl*zl ) 
ff (ind)=fl 
10 continue 


c 

c 

c 

c 


c 


c 

c 

c 

c 


c 


Find worm pitch points 
jt=l 

dr=dabs(rr(l) -r2) 
do 20 j=2,ind 

if (dabs(rr(j ) -r2) .lt.dr) then 
dr=dabs(rr(j ) -r2) 

endif 

20 continue 

if (ff(jt) .gt. endl. or. ff(jt) .lt.startl) goto 25 
fl=ff (j t-1) 

calP root ( f 1 , f r , fnl , yl , aO , ratio , rb , tb ,btr , r2 , epp , f 1 , xl , zl) 
if (fl. gt.endl.and.fi. lt.startl) goto 25 
rain=ds qr t ( ( aO - xl ) * ( aO - xl ) +z l*z 1 ) 


Direction of normal line and relative velocity at pitch point 

call normal ( f 1 , fnl , xl , yl , zl , aO , ratio , rb , tb , btr , 

& px,py,pz,rx,ry,rz) 

write(7 ,6) fnl,xl,zl,px,py,pz 
write (8 , 6) fnl,xl,zl,rx,ry,rz 

nc=nc+l 

tp2(nc)=fn2 

fp2(nc)=fl 

pxj (nc)=xl 

pzj (nc)=zl 

xpn(nc)=px 

ypn(nc)=py 

zpn(nc)=pz 

vx(nc)=rx 

vy(nc)=ry 

vz (nc)=rz 


c 
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25 


do 30 i=l,ind 
xm(i)=0 
zm(i)=0 
rr(i)-0 
ff (i)“0 
30 continue 

50 continue 

return 

end 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


c 


c 

c 


c 

c= 

c 

c 


END of Subroutine PITCH 


Subroutine WORMT called by Subroutine PITCH and ROOT 

This subroutine is used to calculate the coordinates of the worm 
surface in Sj specially for yj=0. 


To run this subroutine the input parameters are: 

aO center distance of the first enveloping (mm) 

btr apex angle of the generating plane (rad) 

ratio speed ratio of the first enveloping 
rb radius of base circle (mm) 

tb tan(bt) 

fnl the given position (tn2*ratio2) (rad) 

ft running parameter of the first enveloping (rad) 


The output of the subroutine are: 

xm.zm the coordinates of the worm profile 

mw the index of the worm profile point on edge curve (mw=l) 


subroutine wo rat ( ft , fnl , y , aO , ratio , rb , tb , btr , xm , zm , mw) 
implicit double precision(a-h,o-z) 


mw=0 

f2=ft/ratio 
fl=f t- fnl 

call eqw(fl,f2,a0,ratio,rb,tb, 

& p , q , r , s , pi , ql , rl , si , p2 , q2 , r2 , s2 ) 

dw=p*(q2*rl-ql*r2)+q*(pl*r2-p2*rl)+r*(p2*ql-pl*q2) 
if(dw.eq.O) goto 10 

fz=- (p*(q2*sl-ql*s2)+q*(pr*s2-p2*sl)+s*(p2*ql-pl*q2) )/dw 
dw=p*rl-pl*r 

xm= (y*(ql*r-q*rl)+(sl*r-s*rl) ) /dw 
zm= (y*(pl*q-p*ql)+(pl*s-p*sl))/dw 
if(zm.lt.fz) goto 10 
goto 20 

10 mw=l 

20 return 
end 


END of Subroutine WORMT 
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c 

c 

c= 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


c= 


Subroutine ROOT called by Subroutine PITCH 

This subroutine is used to search the pitch point from the worm 
profile in Sj specially yj=0. The method used here is called 
golden division. 


To run this subroutine the input parameters are: 

aO center distance of the first enveloping (mm) 

btr apex angle of the generating plane (tad) 

r pitch radius of gear (mm) 

ratio speed ratio of the first enveloping 
rb radius of base circle (mm) 

tb tan(bt) .. 

fnl the given position (tn2*ratio2) (rad) 

fi,fe interval of fl for root searching (rad) 

epp convergence criterion for root searching (rad) 


The output of the subroutine are: 

xl,zl the coordinates of the worm pitch point in Sj 
fl the value of running parameter at pitch point 


subroutine root ( f i , f e , fnl , y , aO , ratio , rb , tb , b tr , r , epp , f 1 , xl , zl ) 
implicit double precision(a-h , o-z) 
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120 

130 

150 

160 

170 


call wo rmt ( f i, fnl, y, aO, ratio, rb ,tb, btr, xm,zm,mw) 

gi=dsqrt( (a0-xm)*(a0-xm)+zm*zm) -r 

call wormt ( f e , fnl , y , aO , ratio , rb , tb , btr , xm , zm , mw) 

ge=dsqrt( (a0-xm)*(a0-xm)+zm*zm) -r 

df-fe-f i 

if (df.lt. epp) goto 150 
f tl=f i+ . 382*df 
ft2=f i+ . 618*df 

call wormt ( f tl , fnl , y , aO , ratio , rb , tb , btr , xm , zm , mw) 
gl=dsqrt( (a0-xm)*(a0-xm)+zm*zm) -r 

call wormt (ft2, fnl, y,a0, ratio, rb, tb, btr, xm,zm,mw) 
g2=dsqrt( (a0-xm)*(a0-xm)+zm*zm) -r 
if (gl*gi.gt.0) goto 120 
fe=ftl 


ge=gl 

goto 110 

if (g2*gi . gt . 0) 

fi=ftl 

fe=ft2 


goto 130 


gi-gl 
ge=g2 
goto 110 
fi=ft2 


g i= S 2 
goto 110 
f 1 = ( f i+f e ) /2 . 0 

call wormt (fl , fnl ,y , aO , ratio , rb , tb ,btr ,xl , zl ,mw) 

goto 170 

mc=l 


return 

end 

— END of Subroutine ROOT 
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c 

c 

c= 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


c 


c 

c= 

c 

c 

c= 

c 

c 

c 

c 

c 

c 

c 


Subroutine NORMAL called by Subroutine PITCH 


This subroutine is used to calculate the propertities of the worm 
pitch point in Sj . 

To run this subroutine the input parameters are: 

aO center distance of the first enveloping (mm) 

btr apex angle of the generating plane ^ (rad) 

fl the value of running parameter at pitch point 

ratio speed ratio of the first enveloping 
rb radius of base circle (mm) 

tb tan(bt) 

xl,zl the coordinates of the worm pitch point in Sj 


The output of the subroutine are: 

wnx,wny,wnz the normal direction components of 
pitch point in Sj 

xw,yw,zw the relative velocity components of 

pitch point in Sj 


subroutine normal ( f 1 , fnl , xl , yl , zl , aO , ratio , rb , tb , btr , 
& wnx,wny ,wnz ,xw,yw,zw) 

implicit double precision(a-h,o-z) 

f2=fl/ratio 
cf2=dcos (f2) 
sf2=dsin(f2) 
c f l=dcos ( f 1 - fnl ) 
s f l=ds in ( f 1 - fnl ) 

wnx=cfl*sf2+tb*sfl 
wny=tb*cf 1 - sf l*sf 2 
wnz=-cf2 

rx=xl*(l-aO/dsqrt(xl*xl+yl*yl) ) 
r y=yl*(l-aO/dsqrt(xl*xl+yl*yl) ) 
rz=zl 

xw=wny*rz - wnz*ry 
yw=wnz*rx - wnx*r z 
z w=wnx*r y - wny *r x 
rlu=dsqrt (xw*xw+yw*yw) 
rld=dsqrt (xw*xw+yw*yw+zw*zw) 
rlead=dacos (rlu/rld) 
return 
end 


END of Subroutine ROOT 


Subroutine LOAD called by MAIN 

This subroutine is used to calculate load sharing among the meshing 
teeth. All the calculations are under the given position tn2. To 
drive this subroutine the subroutines DIMEN, CONTACT and PITCH must 
be called first. 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 


c 


c 

c 

c 


To run this subroutine the input parameters are: 
center distance of wormgear 
the number of real meshing pairs 
number of meshing teeth 
output power from gear shaft 
pitch radius of gear 
speed ratio of wormgear 
input speed from worm shaft 
the position of each meshing pair 
friction coefficient 

the coordinates of pitch points in Sj 
the normal direction components of 
pitch point in Sj 

cont the length of contact lines on each meshing pair 


a2 
nc 

nzmax 

power 

r2 

ratio2 

rpm 

t2 

uf 

xl , zl 

pnx , pny , pnz 


(mm) 


(kw) 

(mm) 

(rpm) 


The other results which are not called by other subroutines are 
written in the data file under the name of 'load.dat' 


subroutine load (a2, ratio2,r2, nzmax, power, rpm, uf,nc,t2, 
& xl , zl , pnx , pny , pnz , cont , wtor , wtorf , f orn , f ornf ) 

implicit double precision(a-h, o-z) 

dimension t2 (nzmax) , cont (nzmax) ,xl (nzmax) , zl (nzmax) , 

& pnx (nzmax) , pny (nzmax) , pnz (nzmax) , 

& rei(10) ,refi(10) ,ww(10) ,wwf (10) 

9 format(2x, 6fl2 .4) 


write (9 ,*) ' 
write(9,*) ' 
write(9 ,*) ' 
&rque ' 

write(9,*) ' 

&- ' 

write(9,*) ' 
&uf O0 ' 
write(9,*) ' 
& ' 


normal load and worm torque distribution' 
pair No . | normal load [ worm | to 


pair No . | uf=0 ufoO uf— 0 


p i=4 . d0*datan ( 1 . dO ) 

grpm=rpm/ratio2 

torqeg=9 . 550*1 . 0e6*power/grpm 

calculation of the sum of the contact lines 

dd=0 . 0 
ddf =0 . 0 
do 20 i=l,nc 

pn=dsqrt (pnx ( i ) *pnx ( i ) +pny ( i ) *pny ( i ) +pnz ( i ) *pnz ( l ) ) 

xn=dabs (pnx ( i ) /pn) 

yn— dabs (pny ( i ) /pn) 

zn=dabs (pnz ( i ) /pn) 

rei(i)=zn*(a2-xl(i)) -xn*zl(i) 

pv=dsqrt(r2*r2+(xl ( i) *ratio2) **2) 

ref i ( i) =uf *r2*r2/pv 

ww(i)=yn 

wwf ( i ) =yn+uf *xl ( i ) *rat io2/pv 
dd=dd+r e i ( i ) *cont ( i ) 
ddf=ddf+(rei(i) -refi(i) )*cont(i) 

20 continue 
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c 

c 

c 


calculate normal load and worm torque distribution (uf=«0) 

do 30 j=l,nc 

gfor=torqeg*cont( j )/dd 
wtor=gfor*ww(j )*xl(j ) 

c calculation of the normal load and worm torque distribution 
c with friction 
c 

gforf=torqeg*cont(j )/ddf 
wtorf=gforf*wwf (j )*xl(j ) 

write (9,9) t2 ( j ) *ratio2 , gf or , gf orf , wtor , wtor f 
30 continue 
c 

return 

end 

c 

c - END of Subroutine LOAD 


130 



Appendix 3 The Principal Relative Curvatures 

Assume that a surface is a simple surface. At a point Q of there 
exists a normal n. The tangent plane of the surface 2^ at Q is perpendicular 
to the normal n. Any plane containing the normal n is perpendicular to the 
tangent plane. This plane is called a normal plane of the surface at point 
Q. The intersection of the normal plane and the tangent plane is a straight 
line. The direction of the straight line is called a tangent direction of 
surface ^ at point Q. Therefore, at a point on there is only one normal 
direction and only one tangent plane, but there are infinitely many tangent 
directions and infinitely many normal planes. 


At point Q the intersection of surface Z x and a normal plane gives a 
planar curve , correspondingly the intersection of the normal plane and the 
tangent plane gives a tangent direction. The curvature of the planar curve at 
point Q is called the normal curvature of surface in this tangent 
direction. Let the point Q be given by the position vector r(u,v), where u, v 
are two independent variables . Then the normal curvature is : 

jc = Ldu 2 +2 Mdudv+Ndv 2 (a.l) 

n Edu 2 + 2 Fdudv+ Gdv 2 


where E=r u 2 , F=r u *r v , G=r v 2 , L=r uu -n, M=r uv -n, N=r w *n; and n is the normal 
direction given by r u xr v / 1 r u xr v | . If the position of Q is given by r(x,y,z) 
where z=z(x,y) represents the equation of the surface, then E, F, G, L, M, N 
can be written as: 


E = 1 + 


, |z )2 

dx 


^ _ dz dz 
dx dy ' 


G = 1 + 


(4p) 2 

dy 


L = 


dPz 
dx? 
DB ' 


M = 


dxdy 
DB ' 


N = 


&z 

dy 2 

DB 


where 


= / 


DB = J 1 + 


<!> 7 




(a. 2) 
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The sign of the normal curvature is determined by the normal direction. 
If the normal direction points to the center of curvature center, the normal 
curvature is positive, otherwise it is negative. 


Along different tangent directions the normal curvatures of a surface 
are different. The principal curvatures are the extreme values of the normal 
curvatures. The corresponding tangent directions are called the principal 
directions (Kreyszig, 1959). 


If two surfaces and S 2 are in meshing at a contact point Q, the 
normals of the two surfaces at the contact point coincide. Along a tangent 
direction let the normal curvatures of the two surfaces be k n <1) and , then 

in this direction the relative curvature of two surfaces is defined as : 


v — It ' 

■‘'•r ~ 


( 2 ) 


(a. 3) 


The principal relative curvatures are the extreme values of the relative 
curvatures , the corresponding directions are the principal relative 
directions. If the two surfaces are simple surfaces, there are two principal 
relative curvatures and two principal relative directions. 


If the two surfaces contact along a line, the direction of the contact 
line is one of the principal relative direction and the corresponding 
principal relative curvature is zero. The other principal relative direction 
is perpendicular to the direction of the contact line and the principal 
relative curvature is (Litvin, 1989; Shen, et al. 1983): 

K = Ei n i~ 2F i m i +g i l i _ e z^2~^F 2 M 2 -*-G 2 L 2 ( a . 4 ) 

E 1 G 1 ~F^ E 2 G 2 - i*2 


In the second enveloping process the gear surface, Eqs. 27, is generated 
by the worm surface, Eqs. 18. In this process the two surfaces are in line 
contact. Let P=A 0 (21) F, the equations of the gear surface can be written as: 
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P = 0 


D^P =0 (a. 5) 

D e P = 0 


By solving Eqs . 18, the relations of z 1 (x 1 ,y 1 ,^(x 1 ,y 1 )) can be found. A point 
on the worm surface can then be expressed as * 1 (x 1 ,y 1 ,z 1 (x 1 ,y 1 )) . Further, by 
solving Eqs. 27, the relations of z 2 (it 2 ,y 2 ) , V>(x 2 ,y 2 ) , 0(x 2 ,y 2 ) can be found. 

A point on the gear surface can be written as r 2 (x 2 ,y 2 ,z 2 (x 2 ,y 2 )) . Then, by 
using Eq. a. 4 the principal relative curvature is obtained. In this equation, 
one should make the following substitutions: 


E, = 


p 2 + p 2 
_ r Xi z i 




F x = 


F v F, 


yi 


a- F y* +F > 


(a. 6) 


L,= 



2 

^ Pz, 


F Vx F Zi 


F n F *> 

2 


M = 



F yi+i F ^i*i 

. A71 = ■ 

^yi*i **iti 


DA 

— / 


DA 

f * ¥ i 

DA 



(a. 7) 


i?, = 


P 2 + P* 

P* 2 

z 2 


P„ P„ 

r* _ *2 ”2 

F 2 


Py 2 2 + P Z2 2 

V 


(a. 8) 


L 2 = m ( £ 2 Pe 2 e 2 “ 2bc P ti e 2+ c 2 P titi ) 

M 2 = m ( bd P e20a -(jbe + cd )P* i02 + ce P +iti ) (a. 9) 

N z = m {cPPqj- 2de P^ i e a + © 2 


where 


DA = ** F, iti 


(a. 10) 
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** 



c = 

** 

P *2 

P *2*l 

P z**i 



p *& 

P *a®a 

p y t 



e = 

p y. 

P *2 

p y*i i 

P «2*l 



P y& 

P *2 Q 2 

p titi 

p <Ma 


m = 

cosP 


P *i®a 

p ®a®a 



gp\ 



(a. 11) 


_ dF F _ &F 

* ~dq ' 97 dgdx 

< q.r = x x , y x , z x , iM 

p _ dP p m &P_ 

P * ds ' 8C ds dt 

( s , t = x x , y x , z x , *1^, x 2 , y 2 , z 2 , 0 2 ) 
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